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ABSTRACT 

In a six-year program started in July 2014, the Extended Baryon Oscillation Spectroscopic Survey 
(eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point 
Observatory. These observations will be conducted simultaneously with the Time Domain Spectro¬ 
scopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA 
Sources (SPIDERS) program designed for studies of X-ray sources. In particular, eBOSS will measure 
with percent-level precision the distance-redshift relation with baryon acoustic oscillations (BAO) in 
the clustering of matter. eBOSS will use four different tracers of the underlying matter density field 
to vastly expand the volume covered by BOSS and map the large-scale-structures over the relatively 
unconstrained redshift range 0.6 < z < 2.2. Using more than 250,000 new, spectroscopically confirmed 
luminous red galaxies at a median redshift z = 0.72, we project that eBOSS will yield measurements 
of the angular diameter distance dA{z) to an accuracy of 1.2% and measurements of H{z) to 2.1% 
when combined with the z > 0.6 sample of BOSS galaxies. With ~ 195, 000 new emission line galaxy 
redshifts, we expect BAO measurements of dA{z) to an accuracy of 3.1% and H{z) to 4.7% at an 
effective redshift of z = 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will 
provide the first BAO distance measurements over the redshift range 0.9 < z < 2.2, with expected 
precision of 2.8% and 4.2% on dA{z) and H{z), respectively. Einally, with 60,000 new quasars and re¬ 
observation of 60,000 BOSS quasars, we will obtain new Lya forest measurements at redshifts z > 2.1; 
these new data will enhance the precision of dA^z) and H{z) at z > 2.1 by a factor of 1.44 relative to 
BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales 
through redshift-space distortion (RSD) measurements, improved tests for non-Gaussianity in the 
primordial density field, and new constraints on the summed mass of all neutrino species. Here, we 
provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral 
quality from early data, and projected cosmological constraints from eBOSS. 

Subject headings: cosmology: observations 
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1. INTRODUCTION 

The origin of the accelerating expansion of the Universe is argnably the most important unknown in physics today 
and has inspired significant efforts to probe beyond the standard model of high-energy physics through observational 
cosmology. The recent measurements of the Cosmic Microwave Background (CMB) from the Planck satellite support 
a picture where the acceleration is driven by “dark energy” with density ^de = 0.692 zb 0.012 in a spatially flat 
universe (Planck Collaboration et al. 2015a). Combining these results with current measurements from baryon acoustic 
oscillations (BAO), Type la SNe, and i^o, the data imply a constant equation of state w = —l-006lo^O9i ^5% 
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confidence, where w is the ratio of pressure to energy density for dark energy. Thus, current observations are generally 
consistent with the simplest picture where dark energy is described completely by Einstein’s cosmological constant 

New precise observations can unravel the origin of the accelerating universe; specifically, to determine if cosmic 
acceleration is caused by deviations in General Relativity (GR) on large scales or by a new form of (dark) energy. It 
is possible to decouple scenarios of acceleration that require dark energy from those that require modifications to GR 
by independently probing both cosmic expansion history and the structure growth rate. Four primary observational 
techniques are generally accepted as the most powerful toward obtaining that goal (e.g. Albrecht et al. 2006): SNe la, 
weak leasing, galaxy clusters, and BAG. Wide-held, optical spectroscopy hgures prominently in three of these probes: 
spectroscopically observed galaxies improve calibration of photometric redshifts for weak leasing; direct spectroscopy 
of cluster galaxies provides precise redshifts and velocity dispersions as a proxy for cluster mass; and spectroscopy of 
galaxies and quasars provides the atlas of large-scale structure in which the BAG feature is embedded. 

The Sloan Digital Sky Survey (SDSS; York et al. 2000) at Apache Point Gbservatory (APG) has consistently provided 
the largest spectroscopic samples for cosmological analysis. In the hrst two generations of SDSS, generally known as 
SDSS-I and -II, redshifts of nearly one million galaxies were measured spectroscopically (Abazajian et al. 2009). The 
Baryon Gscillation Spectroscopic Survey (BGSS; Dawson et al. 2013) performed spectroscopic observations of large- 
scale structure in SDSS-III (Eisenstein et al. 2011); BGSS recently completed spectroscopy on more than 1.5 million 
galaxies as faint as i = 19.9 and more than 150,000 quasars as faint as g = 22. Measurements of BAG with BGSS 
have led to 1-2% precision measurements of the cosmological distance scale for redshifts z < 0.6 and z = 2.5. 

With observations that commenced in July 2014, SDSS-IV will continue this legacy in three distinct surveys. The 
second generation of the APG Galactic Evolution Experiment (APGGEE/APGGEE-2; Majewski et al. 2015) will in¬ 
vestigate the formation and chemical history of the Milky Way using high-resolution, infrared spectroscopy of 300,000 
stars. Mapping Nearby Galaxies at APG (MaNGA; Bundy et al. 2015) will measure the internal structure of approxi¬ 
mately 10,000 galaxies using the BGSS spectrograph with cartridges outfitted for spatially resolved spectroscopy. The 
Extended Baryon Gscillation Spectroscopic Survey (eBGSS) is the new cosmological survey within SDSS-IV. 

The eBGSS program will use the same 1000-fiber optical spectrographs as those in BGSS (Smee et al. 2013). We 
will expand the selection of luminous red galaxies (LRG) beyond that probed by BGSS and obtain better than a 1.0% 
precision distance estimate when combined with the z > 0.6 tail of the BGSS galaxy population. With observations of 
a new sample of emission line galaxies (ELG) over the period 2016-2018, eBGSS will produce a 2.0% precision distance 
estimate at higher redshifts. We will obtain a 1.8% precision distance estimate in the redshift range 0.9 < z < 2.2 using 
quasars that have luminosities and areal densities well-suited to sensitivity of the BGSS spectrographs. Finally, we 
will sharpen the BGSS Lya forest measurements by a factor of 1.44 with a new selection of z > 2.1 quasars, providing 
stronger leverage on the history of dark energy. Goncurrent with the eBGSS survey, we will perform two complementary 
programs that use ~ 100 of the BGSS spectrograph fibers within each field-of-view. The Time Domain Spectroscopic 
Survey (TDSS; Morganson et al. 2015) will target variable stars and quasars. The Spectroscopic Identification of 
eRGSITA Sources (SPIDERS) will target active galactic nuclei, quasars. X-ray emitting stars, and galaxy clusters 
identified in X-ray images of wide-area surveys carried out by RGSAT, XMM-Newton, and, eventually, eRGSITA 
(Merloni et al. 2012). 

With four classes of spectroscopic targets (LRG, ELG, quasar, Lya forest quasar), eBGSS will enable the first high 
precision distance measurements in the epochs when dark energy emerged as the dominant dynamical component 
of the Universe. In addition to BAG distance measurements, eBGSS will provide new tests of GR on cosmological 
scales through redshift-space distortions (RSD), new tests for non-Gaussianity in the primordial density field, and new 
constraints on the summed mass of all neutrino species. 

This paper is one of a series of technical papers describing the eBGSS survey. Details of the LRG target selection 
algorithm is described in Prakash et al. (2015a) while Myers et al. (2015) presents the quasar target selection algorithms. 
The ELG selection is still under investigation; Raichoor et al. (2016), Comparat et al. (2015a) and Delubac et al. 
(2016) present various approaches and results. Here, we provide the scientific motivation behind eBGSS, summarize 
the observation strategy for the cosmological tracers, present the expected quality of spectra for each target class, 
and compute the projected cosmological constraints. In Section 2, we provide a review of the cosmological signature 
of BAG and RSD in the matter power spectrum. We use BGSS measurements to demonstrate the sensitivity of 
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spectroscopic surveys to both signatures and to demonstrate sources of systematic errors that must be addressed in 
eBOSS and future spectroscopic programs. We present the eBOSS program in Section 3, including time allocation, 
expected areal coverage, and requirements for target selection. In Section 4 we discuss the survey strategy, selection of 
the galaxies and quasars that will be used to measure the matter power spectrum, and statistics after fiber assignment. 
Favorable weather during SDSS-III led to an early completion of the BOSS program. A fraction of the remaining time 
was allocated to an eBOSS pilot program known as the Sloan Extended Quasar, ELG, and LRG Survey (SEQUELS). 
Using the results of those 66 dedicated plates, a demonstration of the data quality for each eBOSS target class is 
presented in Section 5. We summarize the expected quality of the data with respect to the low-level specifications 
required to achieve percent-level precision distance estimates in Section 6. We also present plans to improve the data 
reduction software in Section 6. Finally, the cosmological projections for eBOSS are provided in Section 7 and a 
summary is provided in Section 8. 

2. SIGNATURE OF BAO AND RSD IN SPECTROSCOPIC SURVEYS 

Experience from BOSS led to the survey design for eBOSS described in Section 3 and Section 4. BOSS proceeded 
with the primary goal of obtaining new measurements of the cosmic distance scale through BAO at redshifts 2 < 0.6 
and at 2 : = 2.5. While there have been other successful BAO programs (e.g. the 2dE Galaxy Redshift Survey, 6dF 
Galaxy Survey and WiggleZ Dark Energy Survey; Gole et al. 2005; Beutler et al. 2011; Blake et al. 2011b), we only 
review BOSS here because it provides a direct demonstration of spectroscopic constraints on cosmology that we expect 
from eBOSS, and sets the scene for this survey. In what follows, we outline the BOSS observational program in 
Section 2.1, the resulting BAO measurements in Section 2.2, and RSD measurements in Section 2.3. 

2.1. Sample of Galaxies and Quasars from BOSS 

BOSS is described in detail in Dawson et al. (2013), so we provide only a quick summary to highlight the essential 
features that are either inherent or complementary to the eBOSS program. BOSS and eBOSS use the same telescope as 
that used in SDSS-I and -II: the 2.5-meter Sloan Eoundation Telescope at APO in New Mexico (Gunn et al. 2006). The 
BOSS spectrographs (Smee et al. 2013) were built with smaller fibers, new improved detectors, higher throughput, and 
a wider wavelength range than the SDSS spectrographs previously used at APO. There are two spectrographs, each 
covering the wavelength range 361 nm - 1014 nm. The instrument is fed by 1000 optical fibers (500 per spectrograph), 
each subtending 2" diameter on the sky. 

The detection of BAO (Eisenstein et al. 2005) in the SDSS LRG sample (Eisenstein et al. 2001) motivated a similar 
selection of galaxies as the primary sample of spectroscopic targets for BOSS. At the beginning of BOSS observations, 
the techniques for analysis were already advanced to the stage that the collaboration could reasonably expect to perform 
BAO distance measurements at roughly 1% precision with the galaxy sample. In a more experimental program, more 
than 150,000 quasars at redshifts z > 2.1 were selected to measure fluctuations in the matter density field as observed 
through the Lya forest. Contrary to the galaxy program, the Lya forest program was created to make the first 
BAO measurements with a new tracer with an expectation that any detection would be made with somewhat lower 
significance than the galaxy measurement. 

The BOSS survey obtained roughly 10,000 deg^ of spectroscopic coverage over a five year period. At the completion 
of the main BOSS program, fibers were plugged into 2438 unique spectroscopic plates that each cover a circular field 
of view with 3° diameter. All of the BOSS spectra became public in January 2015 with Data Release 12 (DR12; Alam 
et al. 2015a). These data include spectra and classification of 343,160 unique low redshift galaxies (denoted “LOWZ”; 
-^median = 0.32), 862,735 Unique medium redshift galaxies (denoted “CMASS”; 2:median = 0.57), and 181,605 quasars 
between 2.1 < z < 3.5, of which 160,786 are not tagged as broad absorption line (BAL) quasars. 

2.2. Baryon Acoustic Oscillations 

While the overall shape of the power spectrum informs our cosmological model, the most robust cosmological 
measurements from spectroscopic surveys derive from the BAO feature in the clustering of matter. Sound waves that 
propagate in the opaque early universe imprint a characteristic scale in the clustering of matter, providing a “standard 
ruler” whose length can be computed using straightforward physics and cosmological parameters that are constrained 
by CMB observations. The characteristic BAO feature has a comoving scale (roughly 150 Mpc) set by the integrated 
sound speed between the end of inflation and the decoupling of photons around z = 1000. The detection of the acoustic 
oscillation scale (Eisenstein et al. 2005) is one of the signature accomplishments of SDSS. 

Measuring the angle subtended by the characteristic BAO feature at a given redshift provides a means to estimate 
the angular diameter distance, dA^z). Similarly, a measurement of the redshift interval over which the BAO feature 
extends provides a means to directly measure the Hubble parameter H{z) at that redshift. In general, there is a 
combination of dA{z) and H{z) denoted by a generalized distance parameter “R” that is better constrained than 
either of these two quantities alone. In the specific geometry of the measured clustering of BOSS galaxies, the best 
constrained “distance” is approximately R = Dy = [d^iz) cz H~^{z)]^^^ because there are two transverse dimensions 
and one line-of-sight dimension. In the Lya forest, the enhancement of the radial signal gives greater weight to H{z). 

At low redshifts, BAO are a powerful complement to studies with SNe la: they have low systematic uncertainties; 
they measure distances in absolute units whereas supernovae only constrain relative distances; and they can directly 
measure both the cosmic expansion rate H{z) and the distance-redshift relation dA{z). At high redshifts, the large 
comoving volume allows the BAO distance method to obtain remarkably precise measurements of the distance and 
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expansion rate with better precision than SNe la. The BAO method is reviewed in detail by Weinberg et al. (2013, 
see their §4), including discussions of the underlying theory, the effects of non-linear evolution and galaxy bias, survey 
design and statistical errors, control of systematics, recent observational results, and complementarity with other 
probes of cosmic acceleration. 


2.2.1. BOSS Constraints on BAO 

Using data from DR9, the first BAO measurement with the BOSS CMASS sample was published in 2012 (Anderson 
et al. 2012). The DR9 sample covered 1/3 of the final BOSS volume, yet a clear BAO detection allowed a measurement 
of the distance to 2 : = 0.57 with a precision of an = 1.7%. The most recent galaxy clustering measurements with BOSS 
(Anderson et al. 2014b) produce a clear BAO detection using both the LOWZ and CMASS samples of DRll galaxies. 
The BAO detection in the CMASS sample provides a 1.0% measurement of the generalized distance to z = 0.57, the 
most precise distance constraint ever obtained from a galaxy survey. The LOWZ sample produces a measurement to 
2 ; = 0.32 with 2.1% precision. 

A thorough examination of potential systematics in the clustering of the BOSS CMASS galaxies revealed a 10% 
decrease in the detected number density of galaxies when comparing areas with high stellar density to those with 
low stellar density (Ross et al. 2012a). The correlation of galaxy density with stellar density is the most significant 
known bias on measured clustering, likely caused by incomplete deblending of detected objects in crowded fields of 
the SDSS imaging data. On the other hand, no significant correlation is seen between number density and potential 
imaging systematic errors in the LOWZ sample of galaxies. The lack of bias is likely because the lower redshift sample 
appears brighter and is less vulnerable to imaging artifacts such as the deblending problems seen in the fainter, high 
redshift sample. A weighting scheme based on galaxy surface brightness and stellar density was devised to reduce the 
systematic effect of stellar density on the measured clustering signal. When focusing on scales less than 150 Mpc, tests 
on mock catalogs reveal that the weights produce no significant bias on the mean measured correlation function. The 
results imply that systematic errors in galaxy clustering due to artifacts in imaging data (causing density fluctuations 
as high as 10%) can be removed if those artifacts are identified. 

The first measurement of clustering in the Lya forest of BOSS quasars was reported in Slosar et al. (2011). In 2012, 
the first measurements of BAO using the Lya forest detected in the spectra of 2 ; > 2.15 quasars were released. Using 
48,640 quasars in the redshift range 2.15 < 2 ; < 3.5, Busca et al. (2013) detected a peak in the correlation function at a 
separation equal to 1.01 ±0.03 times the distance expected for the BAO peak within a concordance ACDM cosmology. 
Using an alternative analysis of the same quasar sample that included some peculiar quasars and larger wavelength 
coverage, Slosar et al. (2013) reported a 2% measurement error on the distance. In the subsequent DRIl analysis 
(Delubac et al. 2015), 137,562 DRll quasars in the redshift range 2.1 < 2 : < 3.5 were used to detect a peak in the 
correlation function and constrain the distance to a precision of 2.1%. Decomposing the optimized distance, the DRll 
measurement yields a measurement H{z) = 222 ± 7 km s“^Mpc“^ and dA{z) = 1662 ± 96 Mpc at an effective redshift 
0 = 2.34. 

Modeling the cosmic distance scale at 2 : > 2 can be enhanced beyond the Lya forest auto-correlation measurement 
by adding a measurement of the cross-correlation of quasars with the Lya forest absorption field. In particular, the 
relatively high bias of quasars provides leverage to improve precision in the measurement of the angular diameter 
distance relative to the Lya forest auto-correlation. Font-Ribera et al. (2014a) use the DRll sample of quasars to 
measure BAO in the Lya forest-quasar cross correlation both along the line of sight and across the line of sight. 
The best fit correlation function at an effective redshift Ze = 2.36 translates to a measurement H{ze) = 226 ± 8 km 
s“^Mpc“^ and dA{Ze) = 1590 ±60 Mpc. Because the Lya forest sightlines are effectively sampled at random, there is 
no evidence for systematic errors from imaging artifacts in the derived clustering of either Lya forest auto-correlation 
or Lya forest-quasar cross correlation. 

In summary, BOSS provides the most precise BAO measurements to date using luminous galaxies and correlations in 
the Lya forest from high redshift quasars. Highlighting the complementarity between BAO and SNe la, Aubourg et al. 
(2015) explore models with the combination of SNe la and BAO calibrated by the sound horizon scale as measured 
by Planck. They show that the combined Hubble Diagram (Betoule et al. 2014) of SNe la from the Supernova Legacy 
Survey (Conley et al. 2011) and SDSS-H Supernova Survey (Frieman et al. 2008; Sako et al. 2014) can be extrapolated 
to 2 : = 0 when anchored at intermediate redshift by BAO measurements. The BAO distance measurements effectively 
provide calibration of the absolute magnitude of SNe la and lead to a 1.7% precision measurement, Hq = 67.3 ±1.1 
km s“^Mpc“^. This value agrees with the value of Hq derived from CMB data alone under an assumption of flat 
ACDM but is in mild tension with Cepheid-derived distances; for example Riess et al. (2011) report Hq = 73.8 ± 2.4 
km s“^Mpc“^. 

Extending the joint CMB, BAO, and SNe analysis to models of dark energy, Aubourg et al. (2015) report a measure¬ 
ment of the constant equation of state for dark energy w = —0.98 ±0.06 when curvature is treated as a free parameter. 
Allowing an additional free parameter to account for a possible time-evolving equation of state, the joint BAO-SNe 
measurements produce only weak constraints on evolution: Wa = —0.6 ± 0.6, where w{a) = wq + iCa)! — a) (e.g. 
Chevallier & Polarski 2001; Linder 2003). As shown in Figure 1, SNe la provide only weak constraints on the Hubble 
Diagram over the redshift range 1 < 2 : < 2, and BOSS has no sensitivity to the cosmic interval that lies between the 
galaxy and the Lya forest samples. It is in this 0.6 < 2 : < 2 regime where the Universe is expected to transition from 
matter-dominated to dark energy-dominated. New measurements of the cosmic distance scale over this redshift range 
have the potential to improve constraints on models of a time-evolving equation of state for dark energy. 
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Figure 1. Projections for eBOSS LRG, ELG, and quasar distance measurements on a Hubble Diagram presented in comoving distance 
( 77 ) versus redshift. Current BAG measurements from BOSS, SDSS (Xu et al. 2013; Ross et al. 2015), 6 dF Galaxy Survey ( 6 dFGS), and 
WiggleZ (Parkinson et al. 2012) are compared to SNe la measurements (Betoule et al. 2014) and Planck predictions (solid curve) obtained 
by marginalizing over the full likelihood function. 

2.3. Redshift Space Distortions and Modified Gravity 

Dark energy is often invoked to explain current CMB, SNe, and BAO observations that imply an accelerating 
Universe. It is also possible to explain the accelerated expansion of the Universe by modifying gravity at large scales. 
The galaxy redshifts used in spectroscopic BAO measurements of the expansion history help differentiate these two 
possible effects through measurements of the growth of structure via RSD (Kaiser 1987). 

RSD arise because the gravitational pull of matter overdensities causes velocity deviations from the smooth Hubble 
flow expansion of the Universe. These peculiar velocities are imprinted in galaxy redshift surveys in which recessional 
velocity is used as the line-of-sight coordinate for galaxy positions. Although the correlation function of galaxies is 
isotropic in real space, the peculiar velocities lead to an increase in the amplitude of radial clustering relative to 
transverse clustering when the correlation function is measured in redshift space. The resulting anisotropy in the 
clustering of galaxies is correlated with the speed at which structure grows; deviations from GR causing slower or 
faster growth give smaller or larger anisotropic distortions in the observed redshift-space clustering. 

In general, the amplitude of clustering at a given redshift is parameterized by (Ts{z), the rms fluctuations in spheres 
of radius 8 /i“^Mpc. The degree of anisotropy due to RSD depends on the rate of change of the amplitude of clustering. 
This change is typically parameterized as a function of the logarithm of the expansion scale parameter /erg = das/din a, 
where a = (1 + z)~^ is the dimensionless cosmic expansion factor. Because RSD measurements are sensitive to the 
product of the growth rate and the amplitude of matter fluctuations, a wide range in redshift coverage is essential to 
constrain the evolution in clustering amplitude and directly probe gravity. 

2.3.1. BOSS Constraints on RSD and Modified Gravity 

Reid et al. (2012) and Samushia et al. (2013) presented the first measurements and cosmological interpretation of 
RSD in the BOSS DR9 galaxy sample. With these results, they constrain the parameter combination /erg = 0.43±0.07. 
Using the larger DRll sample, Samushia et al. (2014) constrain the parameter combination /tjg = 0.447± 0.028 under 









eBOSS 


7 


an assumption of a spatially flat Universe with cosmological constant. Using the DRIO sample to constrain RSD to 
smaller scales of 0.8-32 h~^ Mpc, Reid et al. (2014) find a model-dependent constraint fas = 0.450 ± 0.011 obtained 
in htting a standard halo occupation distribution model to the anisotropic clustering. Other measurements from the 
same data have found similar results. A summary of current RSD measurements is presented in Figure 2. 

Constraints on gravity from RSD become increasingly powerful as the measurements push to smaller scales, as 
evidenced by the higher precision measurements of Reid et al. (2014) relative to Samushia et al. (2014) and others. 
Cosmological measurements from small-scale clustering are dependent on the accuracy of the modelling on quasi-linear 
and non-linear scales. The development and evaluation of analytic, phenomenological, and halo occupation models for 
anisotropic clustering remains a focus with the BOSS galaxy samples (e.g. Chuang et al. 2013; Beutler et al. 2014b; 
Guo et al. 2015). A study of several models in configuration-space tested against mock galaxy catalogs indicates that 
the clustering signal can be well characterized on scales in the range 40 < s < 80h“^ Mpc (White et al. 2015). Certain 
models, such as those based on Lagrangian perturbation theory, are able to fit the mock clustering samples without 
signihcant bias on scales above 25-30 h~^ Mpc. Continued development of theoretical models that allow use of smaller 
scale data may tighten the current BOSS constraints still further. 
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Figure 2. Current RSD constraints on the growth as a function of redshift compared to the projected measurements from eBOSS. The 
current measurements include those discussed in Section 2.3.1 and those for 6 dFGS (Beutler et al. 2012), the main SDSS sample (Howlett 
et al. 2015), 2dFGRS (Song &; Percival 2009), the SDSS LRG sample (Oka et al. 2014), a recent result from the BOSS CMASS sample 
(Alam et al. 2015b), WiggleZ (Blake et al. 2012), and VIPERS (de la Torre et al. 2013). Various models of modified gravity are shown, 
each with the same background expansion, the same comoving BAO position, and amplitude of the power spectrum normalized to that of 
the CMB at high redshifts. The black curve shows the growth in a ACDM universe, assuming the Planck best fit model parameters. The 
yellow curve shows 7 = 0.5 where / = (Linder 2005). The purple curve shows 7 = 0.6. 



















3. THE EXTENDED BARYON OSCILLATION SPECTROSCOPIC SURVEY 

As described in Section 2, BOSS completed a survey of 10, 000 deg^ and enabled the best BAO distance measurements 
to date at z < 0.6 and the most precise distance measurements of any kind at z ~ 2.5. This success motivated the 
SDSS-IV collaboration to repurpose the BOSS spectrographs for eBOSS. The primary goal of eBOSS is to extend 
the BOSS galaxy measurements to 0.6 < z < 1 and to make the first measurements of clustering with quasars as 
tracers and the first measurements of BAO at redshifts 1 < z < 2. The strategy to make these measurements is not 
purely based on an optimized figure of merit. Instead, we follow a similar philosophy to BOSS and designed eBOSS to 
expand a well-understood sample while exploring new tracers over a fundamentally new redshift regime. We planned 
the expanded galaxy sample to achieve 1% precision on BAO measurements as in the BOSS galaxy sample and the 
new tracers to achieve 2% precision as in the BOSS Lya sample. 

Experience from BOSS led to the design of an LRG sample to measure clustering over redshifts 0.6 < z < 1.0. As 
before, the BAO analysis tools for this type of galaxy are well-established. We designed the eBOSS LRG program to 
match the 1% precision on generalized distance achieved with the BOSS galaxy sample. We use the Fisher matrix 
formalism of Seo & Eisenstein (2007) to determine BAO-based errors on the angular diameter distance and the Hubble 
parameter given an observing volume, number density and bias of galaxies. We assume a monotonically decreasing 
number density and a bias of 1.7 ct8(0)/ct8(z), approximately correct assuming the sample continues the evolution of 
the reddest galaxies observed in BOSS (Guo et al. 2013; Ross et al. 2014). We find that number density can be 
traded for areal coverage for number densities around 50 deg“^ without degrading the precision on the BAO distance 
measurement. This relationship holds true for a fairly wide range of number densities, lending us flexibility in the 
design of the LRG program. Assuming a survey area on the order of thousands of square degrees, the balance of 
density and survey area leads to the general requirement that 300,000 LRG’s in the redshift range 0.6 < z < 1.0 are 
required to achieve 1% precision on the generalized distance. 

Just as BOSS introduced the Lya forest sample as a new probe of clustering, eBOSS will introduce two new selections 
that will expand the legacy of SDSS. The first of these selections will be quasars in the redshift range 0.9 < z < 2.2. 
Because quasars are bright and trace large-scale structure with a high bias, they are the most effective tool to constrain 
BAO at these redshifts. The lower bound on redshift is set to have slight overlap with the eBOSS galaxy samples while 
the upper bound is set to have slight overlap with the BOSS and eBOSS Lya forest sample. The potential of this 
sample is limited by the intrinsic number density of the quasar population. In a series of spectroscopic observations with 
BOSS, 1877 quasars were classified over 14.5 deg^ near the celestial equator. After correcting for targeting efficiency 
and selecting confirmed quasars at 0.9 < z < 2.2 and g < 22, the expected surface density of the quasar population is 
82.6 deg“^ (Palanque-Delabrouille et al. 2013a). We therefore are forced to design the quasar program with a lower 
volume density to sample large-scale structure than the galaxy programs. We make BAO projections assuming a 
flat redshift distribution and the bias relationship b{z) = 0.53 -I- 0.29(1 -I- z)^ determined by Groom et al. (2005) and 
consistent with SDSS measurements (Myers et al. 2007, 2006). Unlike the LRG targets, the limited surface density 
will lead to BAO measurements that are fundamentally limited by shot noise. In order to obtain a 2% measurement 
comparable to the BOSS Lya sample, the quasars must sample the intrinsic population with high completeness over as 
much area as possible. We establish the maximal areal coverage in Section 3.1 and report the required quasar density 
there. 

Potential risk in any one target class can be mitigated by observing a number of different galaxy types overlapping 
in redshift and using cross-correlation between different populations. The predominantly passive galaxies in the LRG 
sample are typically classified by absorption lines against faint continuum emission, leading to a practical limit that each 
target’s z-band magnitude be brighter than roughly zab ~ 20. The emission lines in blue galaxies with significant 
star formation facilitate spectroscopic confirmation at high redshifts with an optical spectrograph. Because ELG’s 
are numerous and have emission line fluxes correlated with the observed u or g-band magnitudes (Comparat et al. 
2015b), they can be selected at a much higher density than the LRG sample. For these reasons, we identified an ELG 
sample over an area that ensures a BAO detection with minimal diversion from the LRG and quasar programs. We 
conservatively assume that the bias of the sample is l.Qa^{Q)/<Js{z), slightly lower than expected (Gomparat et al. 
2013b,a). As with the LRG sample, Fisher projections indicate that the shot noise and areal coverage contribute in 
such a manner that the BAO distance error remains relatively constant for a fixed number of objects for a sample 
with density of hundreds per square degree. We find that a 2% distance BAO measurement can be achieved with a 
total sample size of approximately 190,000 spectroscopically confirmed ELG’s with a uniform redshift distribution in 
the redshift range 0.6 < z < 1.0 or 0.7 < z < 1.1. 

In the DRll BOSS sample, the BAO feature at an effective redshift z = 0.57 was measured at a precision of 0.9% 
using GMASS galaxy targets at a density 85 deg“^ (Anderson et al. 2014a). The BAO feature was measured at an 
effective redshift z = 2.34 at a precision of 2% using Lya quasars at an average density 16 deg“^ (Delubac et al. 2015). 
The marginal value of each Lya quasar therefore exceeds that of a typical galaxy in constraining the matter power 
spectrum at BAO scales, as long as the Lya targets can be reliably identified. We reserved as many fibers as possible 
for efficient detection and observations of Lya quasars for the purpose of improving the Lya forest measurements 
obtained in BOSS. 

Of equal importance to BAO are measurements of RSD with each of these samples. However, as discussed in 
Section 2.3, the current RSD measurement errors are primarily driven by theoretical uncertainties in the modeling 
rather than by the statistical power of the target samples. We therefore rely entirely on the high-level BAO goals 
for determining the observing strategy and properties of the target samples. The requirements for target densities are 
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driven by the goal of obtaining 1%, 2%, and 2% BAO distance measurements on the LRG, quasar, and ELG samples, 
respectively. We also present additional requirements on the uniformity of these samples intended to mitigate potential 
systematic errors in clustering measurements. In this section, we present the program that will meet these high-level 
objectives. 


3.1. The eBOSS Cosmology Program 

Over six years, eBOSS will provide the first percent-level distance measurements with BAO and RSD in the redshift 
range 0.6 < z <2, when cosmic expansion transitioned from deceleration to acceleration. Using LRG, quasars, ELG, 
and Lya absorption as tracers of the underlying density field, eBOSS will probe the largest volume to date of any 
cosmological redshift survey. We designed a program of LRG, ELG, and quasar targets distributed over the available 
fibers with a balance that satisfies the high level goals of obtaining 1-2% distance measurements while also meeting a 
series of practical constraints. 

In BOSS, we typically finished observations of a spectroscopic plate in five 15 minute exposures with 20 minutes 
of overhead due to field acquisition, calibration exposures, and GGD readout. These exposures produced spectra of 
sufficient quality to classify the targets in the GMASS sample as faint as z = 19.9 (AB) with over 98% completeness 
(Bolton et al. 2012). A series of pilot studies using dedicated plates in 2013 and 2014 proved that the same spectral 
depth is adequate to classify quasars from 0.9 < z < 2.2, LRGs to z < 0.8, and ELGs at 0.6 < z < 1.1. Since a typical 
field can be completed in 1.5 hours of exposure and overhead, we anticipate that 1800 plates can be completed over 
the roughly 5400 hours allocated to eBOSS if the 50% weather efficiency experienced during SDSS-III continues for 
SDSS-IV. Considering that each plate will contain 100 fibers dedicated to calibration targets and roughly 100 fibers 
dedicated to SPIDERS or TDSS targets, a total of 1.44M fibers will be available for targets designed to trace large-scale 
structure. Given this expected time allocation and potential sample size, we designed the program as follows: 

1. First, we are constrained by the finite area of a spectroscopic plate (7 deg^) and the requirement of a large-scale 
structure survey to obtain areal coverage with limited gaps between plates. As in previous SDSS programs, each 
eBOSS field center is assigned in the tiling process (Blanton et al. 2003). We chose a density of field centers 
to an average of one plate per 5.0 deg^ to avoid gaps in coverage and mitigate inefficiency in fiber assignments 
caused by density fluctuations in the target sample. At this tiling density, 1800 plates can provide a total survey 
area up to 9000 deg^. The TDSS and SPIDERS projects will each use an average of 50 fibers per plate and 
calibration targets require a minimum of 100 fibers per plate. With 800 remaining fibers per plate, the available 
density of cosmological tracers is 160 deg“^. We use this density of fibers to shape the LRG and quasar samples 
which will be observed together over the majority of the program. We isolated the ELG targets to a series of 
dedicated plates. Those plates will contain almost entirely ELG targets except for a few additional objects at a 
maximum density of 10 fibers deg“^. To conduct the ELG program, we will reduce the maximal area of the LRG, 
quasar, and Lya quasar samples from 9000 deg^ to 7500 deg^, thus leaving 300 plates for the ELG program. 
As discussed below and in Section 4, the density and areal coverage of those ELG targets will be decided based 
on final target selection algorithms and available imaging data. By doing so, we only reduce the volume of the 
primary LRG and quasar samples by 16% but obtain a sample of ELG targets that will provide a BAO constraint 
comparable to the DR9 BOSS galaxy results (Anderson et al. 2012) over a new redshift range. 

2. At an area of 7500 deg^, the Fisher projections predict that 0.9 < z < 2.2 quasars at a density of 58 deg“^ will 
produce a 2% distance measurement. Because the sample is shot-noise limited, any fractional increase in the 
number density translates to the same fractional improvement in the projected BAO uncertainty. Roughly 13 
deg“^ of 0.9 < z < 2.2 quasars are already known from SDSS-I, -II, or -III, leaving only 45 deg“^ quasars to be 
identified and spectroscopically confirmed in eBOSS. Because this sample probes a fundamentally new redshift 
range, we design the program primarily around the observations of quasars, allocating fibers at a density of 90 
deg“^ to ensure that we obtain the minimal target density. All objects assigned to the target class denoted 
QS0_C0RE form the basis of this sample; the purity of this core sample needs to be at least 50% to meet the high- 
level goal of a 2% BAO distance measurement when combined with the sample of known quasars. Throughout, 
purity is defined as the fraction of objects that are reliably classified in the redshift range of interest. 

3. The LRG sample must contain at least 300,000 spectroscopically confirmed galaxies in the redshift range 0.6 < 
z < 1.0 to obtain the high-level goal of a 1% BAO distance measurement. These objects will augment the tail of 
the BOSS GMASS sample in the range 0.6 < z < 0.7 and increase the reach of the luminous galaxy sample to 
z = 1.0. Over 7500 deg^, the corresponding density of galaxies with a successful redshift in the desired redshift 
range is 40 deg“^. In terms of projected BAO errors, a fractional change in the number density is equivalent 
to the same fractional change in area. In BOSS, luminous galaxies were reliably identified in imaging data and 
spectroscopically confirmed in a sample with purity that typically exceeded 95%. The eBOSS sample requires 
higher absolute luminosity to maintain a reasonable rate of spectroscopic classification at the higher redshifts. 
However, the apparent magnitudes are still fainter and features are shifted further into the red, likely reducing 
redshift efficiency. Pilot programs and visual inspection during BOSS originally indicated that we could exceed 
80% purity at 0.6 < z < 1.0. We therefore require that the 40 deg“^ LRG sample be reliably classified and 
at redshifts 0.6 < z < 1.0 from observations of LRG candidates at a density 50 deg“^. As will be shown in 
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Section 4.2 and Section 4.3, the purity of the final sample falls short of requirements, leading to the only failure 
we expect in meeting survey design. 

4. Unlike galaxies, the information extracted from each Lya quasar scales with the S/N of the spectra in the Lya 
forest. Before the beginning of BOSS, McDonald & Eisenstein (2007) presented possible BAO constraints from 
the Lya forest accounting for area, density of quasars, resolution, and S/N. As explained in detail in Section 5.3, 
McQuinn & White (2011) find that the relative value of a quasar roughly doubles as the S/N per A increases from 
one to two, but that deeper spectra approach diminishing returns once the depth exceeds a S/N of three per A. 
To capitalize on the benefits of deeper spectra, we will re-observe z > 2.1 quasars that obtained 0.75 < S/N < 3 
per A in the Lya forest in BOSS. Given that the exposure times in eBOSS will be comparable to those in BOSS, 
we expect the typical S/N to increase by 40% when these quasars are observed a second time. The average 
density of this sample over the BOSS footprint is about 8 deg“^. 

5. After hber allocation to the core quasar, LRG, and repeat Lya forest targets, an average of 12 fibers deg“^ 
remain for additional Lya forest quasar targets. In regions with multiple epochs of SDSS imaging, photometric 
variability was used to enhance the selection of Lya quasars in BOSS (Palanque-Delabrouille et al. 2011). In 
eBOSS, we will use the same technique applied to multi-epoch imaging data from the Palomar Transient Factory 
(PTF; Rau et al. 2009; Law et al. 2009). PTF i?-band imaging data with at least hve epochs covers 90% of the 
available eBOSS footprint. We expect an efficiency of about 20% in the selection of Lya quasars from variability. 
Nominally, we would allocate 12 fibers deg“^ to this sample, but instead allow a target density of 20 deg“^ in 
the PTF regions with an expectation of 3-4 deg“^ new, confirmed Lya forest quasars. We increase the density 
to account for reduced LRG and quasar sample sizes caused by “fiber collisions” between objects that lie too 
close to neighboring targets. 

6. As described above, the ELG program is allocated 300 plates with a goal of obtaining spectroscopic classihcation 
of 190,000 galaxies, where the predicted BAO precision depends primarily on the total number of ELGs rather 
than the density or volume covered. There remains uncertainty in the final ELG program, as the exact selection 
method is still under investigation. The algorithm that preferentially selects galaxies with significant [Oil] A3727 
emission will likely be the one that is most efhcient at obtaining this sample. Given this uncertainty, we require 
only that the sample cover a redshift interval Az = 0.4 to ensure a cosmologically interesting volume, lie at a 
median redshift that is higher than the LRG sample, and have an upper bound 2 ; < 1.1 to avoid confusion in 
[Oil] line identification in the wavelength region that is dominated by sky lines. These constraints effectively 
limit the sample of 190,000 galaxies to lie at either 0.6 < z < 1.0 or 0.7 < z < 1.1. If targets were to HU the 
fibers at the lowest density possible, a sample selected at a density 170 deg“^ over 1500 deg^ would meet the 
requirements if attained at a purity exceeding 74%. Likewise, a sample at 340 deg“^ over 750 deg^ would meet 
the requirements at the same purity. As with the LRG sample, a fractional change in the number density has 
the same consequence for BAO precision as would the same fractional change in areal coverage. 

In total, the final eBOSS spectroscopy will consist of: LRG targets (0.6 < z < 1.0) at a density of 50 deg“^ and 
desired purity exceeding 80%, ELG targets at a somewhat higher redshift over 300 dedicated plates with a desired purity 
exceeding 74%, “clustering” quasar targets (denoted QS0_C0RE) to directly trace large-scale structure (0.9 < z < 2.2) 
at a density of 90 deg~^ and desired purity exceeding 50%, re-observations of faint BOSS Lya quasars (z > 2.1) at a 
density of 8 deg“^, and new Lya quasar candidates identified by variability at an average density of 18 deg“^. 

3.2. Target Selection Requirements 

We will only be able to meet the projected cosmological constraints if our final large-scale clustering measurements 
are limited by statistical errors. This requires stringent control of systematic errors that can modulate the data on 
large scales, such as the impact of stellar contamination and dust extinction on target selection efficiency, variations 
in seeing that alter target selection and redshift success, and so on. These systematics have already been extensively 
studied within BOSS (Ross et al. 2011, 2012a; Ho et al. 2012; Pullen & Hirata 2013), and the greater volume and 
greater statistical power at large scales from eBOSS will place new demands on homogeneity of the target samples. 
These effects are important primarily for measurements (such as neutrino masses, the effective number of neutrino 
species, and inflation) that use the full shape of the galaxy power spectrum. BAO measurements are largely protected 
because they rely on a relatively sharp feature, but with the high precision of eBOSS, one must be careful to extract 
the BAO signature in a way that is immune to large-scale modulations. 

Based on experience from BOSS, we apply a strict list of requirements to the selection of the LRG, ELG, and 
clustering quasar target classes to limit these systematic effects. Because Lya quasars provide random sightlines to 
sample the foreground density field, we do not require uniformity in the Lya quasar sample. 

• Statistical and systematic uncertainties in redshift estimates dampen the BAO feature in the radial direction. 
We require that the data reduction pipeline provide a combined precision and accuracy of (t„ < 300kms“^ RMS 
at all redshifts, where ay is defined as caz/{l + z). Increasing the tolerance in redshift estimates to < 540kms“^ 
RMS would increase the projected BAO and RSD uncertainty by 10%. The identifying features in galaxy spectra 



eBOSS 


11 


are typically well-characterized emission or absorption features with widths less than 300 km s“^. We expect to 
have redshift precision significantly better than this requirement for all classified galaxies; we expect occasional 
line confusion to introduce larger redshift errors that would be classified as “catastrophic”, as addressed in the 
next item. 

The prominent emission lines in quasars are broad and subject to velocity shifts with respect to the systemic 
redshift (e.g. Hewett & Wild 2010; Richards et al. 2011). Measurements of clustering in the BOSS quasar sample 
imply redshift errors as large as A 2 ;/(l + z) = 0.003 (White et al. 2012). In addition to inflated redshift errors, 
Font-Ribera et al. (2013) provides compelling evidence for apparent velocity biases in the 2 : > 2.1 BOSS quasar 
sample of roughly 160kms“^. Comparing automated redshifts to redshifts based on Mgll A2799, we believe we 
can improve the redshifting of quasars from those in BOSS and meet the strict redshift accuracy and precision 
requirements for quasars at z < 1.5. To acknowledge the difficulty in redshift accuracy for the high redshift 
region of the clustering quasar sample, we relax the requirement to [300 -I- 400(z — 1.5)]kms“^ for objects at 
z > 1.5. Any errors from the Mgll redshifts would contribute an additional degradation in the redshift accuracy. 
We estimate that degrading the redshift accuracy by 70% (by the same factor across all redshifts) increases the 
projected BAO and RSD uncertainty in the quasar sample by factor of 1.1. 

• Redshift errors exceeding 1000 kms“^ can potentially bias the derived clustering signal by falsely projecting 
classified spectra onto the wrong redshift, possibly through line confusion or contamination by sky lines. As long 
as the fraction of objects with catastrophic redshift failures is small compared to the number of true tracers at 
any redshift, the misidentified fraction should lead to small enough changes in clustering to not bias the estimate 
of the BAO position. To minimize the impact of catastrophic errors, we require that the spectroscopic pipeline 
classify spectra with fewer than 1% catastrophic errors, where the redshifts are not known to be in error. As 
with the redshift accuracy requirement, we relax the definition of a catastrophic failure in the quasar sample to 
unknown errors in excess of 3000 km s“^. 

• The maximum absolute variation in expected galaxy density as a function of imaging survey sensitivity, stellar 
density, and Galactic extinction must be less than 15% (peak to peak). As discussed in Section 2.2, variations 
smaller than 10% have been mitigated in BOSS by straightforward weighting schemes. Based on this experience, 
we expect that we can extend galaxy weights to regions with 15% variation in target density with little effect on 
measurements in the power spectrum in the range 0.02hMpc“^ < k < 0.3/iMpc“^. Any areas with fluctuations 
greater than 15% could be discarded in the cosmology analysis and thus degrade the final precision of the 
clustering measurement. 

• Finally, for all samples used to directly constrain clustering, we require that our target selection procedures be 
robust against variations in the imaging zeropoints. We require that sample densities vary by less than 15% 
for the estimated uncertainties in zeropoint for the input imaging data in each bandpass. Unlike the previous 
requirement for maximum variation in density, this requirement can be computed directly by evaluating the 
derivative of target density with respect to an offset in photometry in a single band. 

4. SURVEY DESIGN 

Several developments since the beginning of BOSS allow us to identify the new samples of tracers for eBOSS 
cosmological measurements. First, the recent release of infrared sky maps from the Wide-held Infrared Survey Explorer 
(WISE; Wright et al. 2010) makes it possible to expand the selection of LRG targets beyond that probed by BOSS. 
Second, recent spectroscopic observations prove that we can select targets from imaging data and spectroscopically 
conhrm star-forming galaxies with high efficiency from the 2.5-meter Sloan Foundation Telescope (Comparat et al. 
2013b). The galaxies with strong emission lines make it possible to further extend the galaxy redshift survey. Selection 
for these targets is not yet hnalized and we plan to conduct the ELG observations in the third and fourth years of 
eBOSS. Third, spectroscopic observations reveal that quasars in the redshift range 0.9 < z < 2.2 can be efficiently 
identihed from spectroscopy (Palanque-Delabrouille et al. 2013a). Finally, improved photometric target selection 
and time-domain imaging data provide a dense sample of z > 2.1 quasars to signihcantly sharpen the cosmology 
measurements from the Lya forest relative to BOSS. 

The procedures to select targets and complete a program of spectroscopy are well established from BOSS. We will 
filter the photometric catalogs, design the spectroscopic plates, and acquire the spectroscopic data in a nearly identical 
fashion to the procedures outlined in Dawson et al. (2013). However, there are changes to each step in the process, 
listed below. 


4.1. Imaging and Catalog Creation 

The baseline data set for eBOSS target selection is the well-understood photometry obtained from the SDSS imaging 
camera (Gunn et al. 1998) in ugriz filters (Fukugita et al. 1996). The median 5-a depth for photometric observations 
of point sources is m = 22.15, g = 23.13, r = 22.70, i = 22.20, z = 20.71. Photometry for each filter X is presented in 
PSF magnitudes (denoted Apsf), fiber2 magnitudes (Afib 2 ) to represent the Hber aperture losses, model magnitudes 
(Aniodei), or cmodel magnitudes (Acmodei), as discussed in online SDSS documentation^. The SDSS photon counts 

^ http://www.sdss.org/drl2/algorithms/magnitudes 
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and object detections used in eBOSS target selection algorithms stem from the DR9 public release (Ahn et al. 2012). 
The BOSS DR9 photometry is internally calibrated using the “ubercalibration” procedure described in Padmanabhan 
et al. (2008), with residual systematic errors of approximately 1.5% in u and 1% in the other four bands. eBOSS has 
recalibrated the flat fields in all five filters and zero-points in the griz bands using a combination of PanSTARRS-1 
(Kaiser et al. 2010) and SDSS stellar photometry. The residual systematics are reduced to 0.9, 0.7, 0.7 and 0.8% in 
the griz bands, respectively (Finkbeiner et al. 2015). In addition, some poorly-constrained zero-points with errors 
exceeding 3% in the DR9 data are now significantly improved. 

The photometric selection for the BOSS galaxy sample was based entirely on imaging from SDSS. Photometry from 
the UKIRT Infrared Deep Sky Survey (UKIDSS; Lawrence et al. 2007) and the Galaxy Evolution Explorer (GALEX; 
Martin et al. 2005) were used to enhance the selection of Lya forest quasars (Bovy et al. 2012). The use of additional 
imaging resources to identify the primary cosmological sample set a precedent that we follow in eBOSS. The primary 
source of external imaging data will be the highly-uniform, infrared (IR) photometry from the WISE satellite. The 
catalog of IR sources is derived using forced photometry informed by positions of sources in the SDSS imaging data. 
The full process is explained in Lang et al. (2014). The resulting 3.4/im and 4.6/rm {W1 and 1F2, respectively) 
magnitudes can be found in the public release^. 

The imaging sample to be used for selection of ELG targets is still under investigation. ELG targets could be selected 
using SDSS griz and {/-band imaging from the South Galactic Gap U-band Sky Survey (SGUSS) conducted at the 
2.3m Bok Telescope on Kitt Peak. The reductions of those two imaging samples are complete and photometric catalogs 
are available over the entire SGG area. ELG targets could also be selected from deeper grz imaging obtained from the 
Dark Energy Camera (DECam; Flaugher et al. 2012). The deeper imaging data allows more precise identification of 
strong emission-line galaxies in the redshift range of interest. We currently have catalogs based on preliminary data 
reductions of first year DECam observations. Through a large program approved for the Dark Energy Spectroscopic 
Instrument (DESI; Levi et al. 2013), we expect to cover the majority of the SGG with DECam in grz. Known as the 
DECam Legacy Survey (DECaLS^), the raw data from those observations are immediately public and the DESI team 
plans to publicly release full reductions of these data. However, observations are scheduled over 2014 - 2017 and it 
is not yet clear how much area will be available for selection of eBOSS targets on the timescales needed for SDSS-IV 
observations. 

Einally, we use multi-epoch imaging data from the Palomar Transient Eactory (PTF; Law et al. 2009) and data from 
the Faint Images of the Radio Sky at Twenty-Centimeters (FIRST; Becker et al. 1995) program to increase the number 
of quasars. PTF covers roughly 90% of the potential eBOSS footprint at a depth that is sufficient to produce and 
characterize lightcurves. These data are made publicly available by the PTF team^. The average field is observed over 
four years in 45 60-second exposures. Lightcurves are constructed with a customized image processing pipeline that 
produces coadditions of all images acquired in a single year. The depth of a resulting epoch is typically comparable to 
the depth of SDSS photometry. Earlier SDSS photometry is used to anchor the lightcurves and provide a longer time 
baseline to measure variability. The catalog of sources from EIRST is used to identify potential quasars matched to 
SDSS optical counterparts. Sources are selected from the June 5, 2013 version of the FIRST point source catalog®. 

For all imaging samples, magnitudes are corrected for Galactic extinction using the Schlegel et al. (1998) models of 
dust absorption. Galactic extinction coefficients have been updated as recommended in Schlafly & Finkbeiner (2011). 
The extinction coefficients i?„, i?g, Rr, Ri, and Rz are changed from the values used in BOSS (5.155, 3.793, 2.751, 
2.086, 1.479) to (4.239, 3.303, 2.285, 1.698, 1.263), respectively. An extinction coefficient Rwi = 0.184 is used for the 
WISE 3.4/rm band and Rw 2 = 0.113 is used for the 4.6/xm band (Eitzpatrick 1999). 

In BOSS, we used a series of imaging flags to remove possible contaminants from the target samples. The standard 
GALIB_STATUS SDSS imaging flags (SATUR etc.) become increasingly poor at distinguishing real objects from artifacts 
at fainter magnitudes. We therefore make no explicit cuts on our galaxy or quasar targets based on the standard 
SDSS imaging flags. We do remove quasar targets with IMAGE_STATUS flags identified as “BAD_ROTATOR”, 
“BAD_ASTROM”,“BAD_FOGUS’, “SHUTTERS”, “EFJETALS”, “DEAD_GGD”, or “NOISY_GCD”. We rely on our 
use of multi-wavelength photometry (such as imaging from WISE) to identify true astrophysical objects as candidate 
sources for spectroscopy. We expect that some objects close to bright stars might have contaminated colors, causing 
spurious objects to enter into the target selection. Areas close to bright stars represent a tiny fraction of the sky and 
it is not obvious a priori how close to bright stars the photometry remains reliable, so we do not explicitly remove 
targets. Glustering analyses will likely either mask or re-weight the galaxy and quasar densities due to the loss of 
objects near bright stars, as has been done for the BOSS clustering analyses (e.g. Anderson et al. 2012). 

4.2. Target Selection Algorithms 

The targets for eBOSS spectroscopy will consist of: LRGs at 0.6 < 2 < 1.0, ELGs at somewhat higher redshift, 
“clustering” quasars to directly trace large-scale structure (0.9 < z < 2.2), re-observations of faint BOSS Lya quasars 
(z > 2.1), and new Lya quasars (z > 2.1). The selection algorithms for identifying the LRG and quasar samples are 
complete and observations of these targets have begun. The selection algorithm for ELG sample is under investigation. 
Here, we summarize the selection scheme and statistics of each sample. 

^ http://unwise.me/ 

^ http://legacysurvey.org/ 
http://irsa.ipac.caltech.edu/Missions/ptf.html 

® http://sundog.stsci.edu/first/catalogs/readme_13jun05.html 
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Table 1 

Expected redshift distribution for the LRG, quasar, and Lyo quasar samples. The surface densities are presented in units of deg“^ 
assuming that 100% of the objects in the parent sample are spectroscopically observed. Entries highlighted in bold font denote the 
fraction of the sample that satisfies the high-level requirement for the redshift distribution of the sample. 



GMASS 

LRG 

LRG 

qS0_C0RE 

QSQ_CQRE 

Lyo: Quasar 

Quasar 


Known 

(Zconf > 0)"" 

(■^conf ^1) 

New 

Known'^ 

Reobservation 

PTF'^ 

Poor Spectra 

- 

4.0 

6.7 

3.2 

0.9 

0 

0.7 

Stellar 

- 

5.3 

5.3 

2.8 

0.3 

0 

10.7 

Galaxy 

- 

- 

- 

6.6 

0.4 

0 

1.5 

0.0 < z < 0.5 

27.3 

0.6 

0.6 

1.0 

0.4 

0 

0.2 

0.5 < 2 < 0.6 

45.7 

6.2 

5.9 

1.1 

0.4 

0 

0.1 

0.6 < z < 0.7 

19.4 

15.2 

14.8 

1.4 

0.7 

0 

0.2 

0.7 < z < 0.8 

3.5 

15.3 

14.7 

1.4 

1.3 

0 

0.2 

0.8 < z < 0.9 

0.2 

9.4 

8.7 

2.2 

1.5 

0 

0.2 

0.9 < z < 1.0 

0.03 

3.2 

2.7 

3.6 

1.0 

0 

0.3 

1.0 < z < 1.2 

0 

0.6 

0.5 

8.4 

1.8 

0 

0.5 

1.2 < z < 1.4 

0 

0 

0 

10.3 

1.8 

0 

0.4 

1.4 < z < 1.6 

0 

0 

0 

10.3 

2.1 

0 

0.6 

1.6 < z < 1.8 

0 

0 

0 

9.9 

2.0 

0 

0.4 

1.8 < z < 2.0 

0 

0 

0 

9.2 

1.9 

0 

0.3 

2.0 < z < 2.1 

0 

0 

0 

4.0 

1.0 

0 

0.2 

2.1 < z < 2.2 

0 

0 

0 

2.2^ 

1.6'^ 

0.5 

0.2 

2.2 < z < 2.4 

0 

0 

0 

1.8'^ 

4.5'^ 

2.9 

0.5 

2.4 < z < 2.6 

0 

0 

0 

1.1^ 

3.1'^ 

1.9 

0.5 

2.6 < z < 2.8 

0 

0 

0 

0.7^ 

1.4^ 

1.0 

0.6 

2.8 < z < 3.0 

0 

0 

0 

0.3'^ 

p 

bo 

0. 

0.7 

0.5 

3.0 < z < 3.5 

0 

0 

0 

0.4^ 

1.2^ 

0.9 

0.5 

z > 3.5 

0 

0 

0 

0.1'^ 

0.1'^ 

0.3 

0.2 

Total Targets 

23 

60 

60 

81.8 

30.4 

8.3 

20 

Total Tracers 

23.1 

43.1 

41.0 

57.9 

13.1 

8.3 

3.2 


^The redshift distribution for the LRG sample is determined by visual inspection. It is not clear how well the automated pipeline will 
perform in the final analysis, so we include the results from visual inspections with ^conf > 0 in the first entry as the more optimistic 
estimate and results with Zconf > 1 in the second entry as the less optimistic estimate. 

'^Objects that satisfy the QS0_C0RE selection that have reliable spectra from previous incarnations of SDSS will not be assigned fibers in 
eBOSS. 

'^The density of PTE-selected quasars is only computed over areas where targets are selected. Roughly 10% of the eBOSS footprint will 
not have PTF-selected quasar targets. 

‘^The population of z > 2.1 quasars that appear in the target selection for the clustering sample will be used for Lya forest studies. 

4.2.1. LRG Samples 

A full investigation of the LRG selection from imaging data is presented in the companion paper (Prakash et al. 
2015a). The final sample was designed to extend the BOSS galaxy sample to higher redshifts following initial studies 
of an optical and infrared selection of LRGs (Prakash et al. 2015b). LRG candidates will be observed (assigned a fiber) 
at an average density of 50 deg“^ as governed by the survey design explained in Section 3.1. Targets are selected at a 
density to oversubscribe the fiber budget and ensure a high overall efficiency of fibers assigned to cosmological tracers. 
Given the fiber assignment statistics described in Section 4.3, we find that a parent sample with density of 60 deg“^ 
is adequate to satisfy the requirement to spectroscopically observe objects at a density of 50 deg“^. 

At the redshifts of the LRG sample {z > 0.6), the 4000 A break is shifted into the SDSS i filter, increasing the 
error of ugr photometric estimates. In addition, the degradation of photometry at higher redshifts provides motivation 
to include ancillary photometric information beyond SDSS to improve the selection efficiency. The WISE Wl filter 
centered at 3.4 microns provides restframe coverage of the “1.6 micron bump” that results from a reduction in the 
opacity of H“ ions exhibited by old stellar populations (John 1988). The final LRG selection algorithm is tuned to 
identifying objects that are red in r — i, i — z and r — Wl colors as demonstrated in Prakash et al. (2015a). 

As will be explained in Section 5.1, the estimate of the redshift distribution stems from visual inspections that are 
difficult to translate to the final pipeline performance. For that reason, we report maximum and minimum number 
densities that correspond to more or less optimistic interpretations of the visual inspections. Between 68.3% and 71.8% 
of the targets identified by this selection lie between 0.6 < z < 1.0, depending on which level of confidence is assumed 
in the visual inspection process. The surface density of the sample is shown as a function of redshift in Table 1. 
Objects with spectra that fail to produce a reliable classification are labeled “Poor Spectra”. 

The density of M-stars and of galaxies at 0.5 < z < 0.6 are both higher than what was expected following the initial 
pilot programs in BOSS. Because of these two contaminants, the parent population represented by the SDSS/WISE 
selection fails to meet the goal of 80% purity for galaxies in the redshift range 0.6 < z < 1.0. Roughly 24% of the 
GMASS galaxies in BOSS lie at redshifts z > 0.6 and roughly 4% lie at redshifts z > 0.7. While not originally intended, 
those GMASS galaxies can be combined with the new eBOSS LRG galaxies to increase the overall sample size and 
obtain a 1% precision distance measurement. 

We assess uniformity of the target sample through a regression analysis of surface density against tracers of potential 
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systematics. We focus on systematics associated with imaging data and astrophysical effects such as dust extinction 
and stellar density. The resulting regression fits reveal that 92% of the potential imaging area in eBOSS has predicted 
surface density that varies by less than 15%, thus satisfying the imaging uniformity requirements outlined in Section 3.2. 
The 8% of the area that fails to meet this requirement must be re-examined when spectroscopic observations are 
complete to assess variations after identification of contaminating sources such as stars. Likewise, we assess the 
variations in number density associated with errors in the zeropoint calibration. We find that zeropoint errors of 0.01 
magnitude in the rizWl filters cause fractional changes in the number density of LRG targets of 2.26%, 2.5%, 6.24%, 
and 0.6%, respectively. Likely due to the fact that it is used for both a color boundary and a magnitude boundary, the 
target selection is most sensitive to the uncertainty in the z-hand calibration; non-uniformity with 15% peak-to-peak 
amplitude occurs in regions where the z-hand zeropoint is in error by ±0.012 magnitudes. Assuming a Gaussian 
dispersion of 0.008 magnitudes in the z-band calibration as discussed in Section 4.1 and in Finkbeiner et al. (2015), 
13.3% of the eBOSS volume exceeds the 15% uniformity requirement. This effect is not expected to be correlated 
with imaging systematics and in the worst case scenario, will have to be addressed through an independent census of 
zeropoints in each SDSS field. Analysis through mock catalogs will inform the manner in which zeropoint uncertainties 
are modeled in the clustering measurements. The full methodology of the uniformity tests is presented in Prakash 
et al. (2015a). 


4.2.2. Quasar Samples 

The algorithm to identify quasars in BOSS was intended to build a sample of z > 2.1 quasars to map the large 
scale distribution of neutral hydrogen via absorption in the Lya forest. As explained in Ross et al. (2012b), the most 
uniform quasar sample for BOSS was identified from a density of 20 targets deg“^ using the “Extreme Deconvolution” 
(XDQSO; Bovy et al. 2011a,b) selection. This selection was performed after estimating the relative density of stars 
and quasars as a function of color, magnitude, and photometric uncertainty. The probability that an object is a quasar 
is determined by the fraction of objects with similar photometric properties expected to be a quasar and not a star. 

Following that precedent, the XDQSO algorithm will be again used to identify quasars for the QSD_C0RE sample in 
eBOSS, but using the XDQSOz version of the algorithm (Bovy et al. 2012) that can be applied to select quasars in 
any redshift range. 

To improve the completeness of the sample, objects with much lower XDQSO probabilities will be included in eBOSS 
relative to BOSS. The inclusion of lower probability objects will also tap into the population of z > 2.1 quasars that 
were missed in BOSS and enhance the sample for Lya forest clustering measurements. 

Colors between optical and WISE passbands can be used to distinguish quasars from stellar objects characterized by 
blackbody spectra (for example Stern et al. 2012). Morphology selection helps reduce the number of lower luminosity, 
extended AGN sources in favor of compact, point-like quasars which map onto higher mass halos. In these ways, WISE 
imaging will be used to decrease stellar contamination and increase the overall efficiency of fiber assignments. 

The final selection algorithm for quasars in the QS0_C0RE sample is described in the companion paper on quasar 
target selection (Myers et al. 2015). The hnal selection when applying XDQSOz and the SDSS/WISE colors results 
in a target density of 115 deg“^. The average density of targets falls to roughly 90 deg“^ after removing objects that 
were confidently classified as a star, galaxy, or quasar in SDSS or BOSS spectra. 

As shown in Table 1, pilot studies in BOSS indicate that the selection identifies 71 quasars deg“^ over 0.9 < z < 2.2, 
of which 13 will already be known from previous observations. The selection will also increase the number of Lya 
quasars: 6.6 new z > 2.1 quasars deg“^ are expected. In addition to stars, we list the rate of galaxy contamination in 
the quasar sample; roughly 8% of the new targets turn out to be galaxies with significant line emission. If targeted at 
100% completeness, this shot noise-limited sample would exceed the goal of 58 deg“^, leading to potential improvements 
of up to 20% over the high-level goal of a 2% precision distance measurement with this sample. 

Performing the same regression analysis as on the LRG sample, we find that ~ 90% of the eBOSS area satisfies the 
requirement of < 15% peak-to-peak variation in the quasar target density. The quasar target selection is very robust 
against zeropoint errors. We find that zeropoint errors of 0.01 magnitudes in any of the ugrizWlW2 filters cause 
fluctuations in target density of less than 1%. The largest contributor to density fluctuations is the (/-band calibration; 
a 0.01 magnitude error leads to a 0.86% change in target density. Given the expected RMS error in flux calibration is 
only 0.009 magnitudes, we expect an RMS scatter in target density of only 0.77% due to calibration errors. Only a 
negligible fraction of sky will see density fluctuations as large as 15%. 

Quasars selected for Lya forest studies are not subject to the same strict requirements of uniformity as those used in 
direct clustering work. Eor this reason, a series of selections was used to increase the total number of z > 2.1 quasars. 
The first selection was based on known quasars from BOSS. A known quasar is included in the eBOSS selection if it 
obtained a low signal-to-noise in the BOSS observation. The mean signal-to-noise ratio is computed over the range 
1040 < A < 1200 A. Objects that did not have broad absorption lines identified in visual inspection and that have 
0.75 < S/N < 3 or S/N = 0 are observed again in eBOSS. As explained in Myers et al. (2015) and in Table 1, we 
also identify unique new Lya forest quasars in the PTF data at a density of 3.2 deg“^ where PTF imaging is available 
to generate sufficient lightcurves. At redshifts z > 2.5, the PTF selection identifies quasars at a density 2.7 deg“^, of 
which 2.3 deg“^ are unique to variability selection. Finally, a small number of objects from the FIRST catalogs are 
included as possible quasars. Because the target density is only around 1 deg“^, we do not include the statistics from 
that sample in Table 1. 
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4.2.3. ELG sample 

The ELG survey will begin in Fall 2016, the third year of eBOSS observations. By starting the ELG program two 
years later than the LRG and quasar programs, the eBOSS team has given itself time to perform deeper analysis of 
the potential selection algorithms. We conducted a series of observations in Fall 2014 to test possible techniques for 
the selection of ELG targets. At the time of writing, the tests are being used to evaluate the redshift success rates, 
redshift distributions, and rates of stellar contamination from four different selection algorithms. 

Each of these four options relies on selecting targets from a subset of imaging data selected from SGUSS [/-band, 
SDSS griz, WISE Wl and W2, or grz imaging with DEGam. The tests will also inform the final choice for imaging 
data to be used in the ELG selection. The expected redshift distribution for each candidate selection is presented in 
Table 2. 


Table 2 

Expected redshift distribution for the four potential ELG selection algorithms. The surface densities are presented in units of deg“^ 
assuming that 100% of the objects in the parent sample are spectroscopically observed. Entries highlighted in bold font denote the 
fraction of the sample that satisfies the high-level requirement for the redshift distribution of the sample. 



gri -\- Uri 

Fisher Discriminant 

DEGam® (Low Density) 

DEGam® (High Density) 

Poor Spectra 

64.7 

19.8 

17.1 

24.1 

Stellar 

4.0 

1.5 

0.0 

0.0 

0.0 < 2 ; < 0.1 

2.1 

2.9 

0.6 

0.6 

0.1 < 2 : < 0.2 

2.7 

2.0 

1.7 

2.4 

0.2 < 2 : < 0.3 

3.6 

2.5 

1.6 

2.0 

0.3 < 2 : < 0.4 

4.1 

1.7 

1.0 

1.2 

0.4 < ^ < 0.5 

4.8 

2.8 

0.6 

0.7 

0.5 < ^ < 0.6 

9.4 

7.2 

1.2 

1.7 

0.6 < ^ < 0.7 

27.6 

25.4 

3.3 

3.7 

0.7 < ^ < 0.8 

42.8 

46.2 

40.7 

44.1 

0.8 < 2 ; < 0.9 

25.9 

40.4 

65.3 

74.1 

0.9 < 2 : < 1.0 

10.8 

17.0 

30.9 

43.7 

1.0 < ^ < 1.1 

2.5 

6.5 

11.1 

16.9 

1.1 < ^ < 1.2 

0.4 

2.6 

5.1 

7.7 

2 > 1.2 

2.2 

2.9 

9.5 

16.9 

Total Targets 

203.9 

182 

190.4 

241.3 

Total Tracers 

107.1 

129.0 

148.0 

178.8 


^The desired range for DECam-selected ELGs is 0.7 < 2 : < 1.1 instead of 0.6 < 5; < 1.0. 


The first candidate for ELG selection uses SGUSS imaging for deeper [/-band photometry to complement SDSS gri 
photometry. Objects for the spectroscopic sample are taken from the union of a g — r and r — i selection in SDSS 
imaging and aU — r and r — i selection in SGUSS and SDSS imaging. In addition to several cuts designed to reject point 
sources and ensure good photometry, objects satisfying g — r < 0.8 and r — i > 0.8 are included in the SDSS gri color 
selection. These cuts lead to a mean target density of 132.5 deg“^. Objects satisfying i — U> —3.5 * {r — i) + 0.7 and 
r — i > 0.7 are included in the SGUSS and SDSS Uri selection at a density of 84.0 deg“^. In all cases, photometry is 
assumed in extinction-corrected model magnitudes. The selections are described in detail in Comparat et al. (2015a). 
The combined algorithm allows the selection of 180 objects per square degree over any region of the South Galactic 
Gap. The density in the test region was somewhat higher and had median redshift 2 ; = 0.72. Many objects in this 
selection lie near the SDSS detection limit; photometric errors are sometimes as large as 0.5 magnitudes. The large 
errors reduce the precision of the color estimates and effectively broaden the redshift distribution. Only 52.5% of the 
targets produce reliable redshift estimates in the redshift range intended for clustering studies. In addition, the large 
photometric errors allow fainter objects into the selection, thus reducing the bias relative to a sample represented by 
a more precisely defined magnitude limit. 

The second candidate selection relies on the addition of the WISE photometry to the SDSS and SGUSS photometry. 
This scheme imposes a cut on Fisher discriminant quantities (Fisher 1936) instead of cuts in color-color diagrams. In 
this case, the Fisher discriminant quantity is a linear combination of colors taken from UgrizWl. The exact definitions 
of the linear combination of colors, the training from a spectroscopic sample, and the expected results are described in 
Raichoor et al. (2016). The selection algorithm can be used to identify targets up to a density of 200 deg“^, somewhat 
higher than the density in the test region. The test observations demonstrate a significant improvement over the 
SDSS-I-SGUSS {gri + Uri) selection. Roughly 71% of these objects are expected to produce reliable redshifts in the 
range 0.6 < 2 : < 1.0. The initial tests of this selection approach the requirement of 74% purity presented in Section 3.1. 

The final two candidates for ELG selection use DEGam grz photometry instead of SDSS, SGUSS, or WISE photome¬ 
try. Both of these selections exceed the 74% purity requirement as detailed in Gomparat et al. (2015a). The advantage 
of these selection schemes is that deeper DEGam photometry allows more precise selection of star-forming galaxies at 
somewhat higher redshifts than the SDSS and SGUSS selections. For this reason, we increase the required redshift 
range for the spectroscopically-confirmed tracer population from 0.6 < 2 ; < 1.0 to 0.7 < 2 : < 1.1. The disadvantage of 
this selection scheme is that it relies on a relatively new imaging camera and a reduction pipeline that is in a state of 
development. The two candidate algorithms were tuned from spectroscopic observations of a broad color selection that 
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produced targets to g < 22.8 at a density of 597 deg“^ in the test region. The DES imaging is reduced using magnitudes 
determined in a 4" diameter aperture and in photometry similar to model magnitudes from SDSS. Extinction-corrected 
magnitudes are used throughout. Both selection algorithms apply the following cuts to reduce stellar contamination 
to nearly zero: gAPER,4 — Sdetmodel < 2, tapera — ?'detmodel < 2, and zaper,4 — ^detmodel < 2, where dfAPER,4 
is derived from an aperture two arcseconds in diameter and -^detmodel is a DES measurement determined from 
SExtractor (Bertin & Arnouts 1996) similar to modelmag in SDSS. The point source rejection also limits the number 
of quasars entering into the sample; the expected quasar contamination is expected to be less than 1% at z < 1.1. 

The first of the DECam selections (denoted “Low Density”) was tuned to achieve a target density of roughly 190 
deg“^. An equivalent to model magnitudes is used in the selection. The color cuts that define that sample are as 
follows: 

• 21.6 < g < 22.8 

• 0.3 < g — r < 0.7 

• 0.25 < r — z < 1.4 

• r — z > 0.45 * {g — r) + 0.4 

• r — z < 0.7 * {g — r) + 0.8 

The resulting redshift distribution is shown in the fourth column of Table 2. 77.8±1.1% of objects lie at 0.7 < z < 1.1 
with a median redshift z = 0.86. 

The second of the DECam selections was tuned to achieve a target density of roughly 240 deg“^. The color cuts 
that define that sample are as follows: 

• 21.5 <g< 22.8 

• 0.2 <5-r < 0.7 

• 0.25 < r — z < 1.4 

• r — z > 0.45 * {g — r) + 0.4 

• r — z < 0.8 * (g — r) -I- 1 

The redshift distribution from this “High Density” DECam selection is shown in the last column of Table 2. 74.1±0.5% 
of objects lie at 0.7 < z < 1.1 with a median redshift z = 0.87. 

While it is possible to increase the density above 240 deg“^ using DECam data, initial tests reveal that contamination 
from galaxies outside the desired redshift range reduces the efficiency of the selection. In addition, as the selection 
approaches 300 deg“^, the density peaks around 3 x 10“^ h^ Mpc“^. At this density, BAO measurements become 
dominated by sample variance rather than shot noise and observing time is more efficiently spent expanding the survey 
volume than by increasing the density. When comparing the low and high density DECam selections presented above, 
one must also consider the strength of the [Oil] and [OIII] emission line fluxes. The weighted mean of the [OIII] line 
flux is 6.6 and 6.9 x 10“^^ erg cm^ s“^ for the low and high density selections, respectively. Likewise, the [Oil] line 
flux is 8.15 and 8.5 x 10“^^ erg cm^ s“^. Although the high density selection produces a somewhat smaller rate of 
galaxies at 0.7 < z < 1.1, the typical line strengths are 5% higher than in the low density selection, thus making it 
more robust to automated classification. 

Tests of uniformity, sensitivity to zeropoint uncertainty, and average target density are underway (Delubac et al. 
2016). A final decision on the selection to be used for eBOSS is expected in early 2016. The sample that is able to 
produce a uniform target density, redshift classification exceeding 74% efficiency over the appropriate redshift range, 
and imaging area sufficient for a total sample of 300 plates will be chosen. If more than one selection meets these 
requirements, the selection algorithms producing the highest median redshift will be the one used for the ELG sample. 

4.3. Tiling and Fiber Assignment 

The goal of survey tiling is to create a spatial distribution of tiles that maximizes the number of targets observed 
with the minimum number of tiles. We define the tiling completeness as the fraction of objects in a given class that 
were assigned fibers. We refer to the fiber efficiency as the fraction of available fibers used for unique science targets. 
In BOSS, the quantity of highest priority was the tiling completeness for ‘decollided’ targets. The physical size of the 
ferrules that support each fiber in the plug plate limits the proximity of neighboring targets to 62". Groups of targets 
that lie within 62" of one another are denoted “collision groups”. The decollided set contains all targets that are not 
within collision groups, combined with the subset of collided targets that can be assigned fibers on a single plate. A 
collision pair contributes one galaxy to the decollided set because, in all cases, one target from the pair will be assigned 
a fiber. If the pair is located within a region observed by more than one tile, the second object may be assigned a fiber 
as well. Thus the completeness of the collided set will be non-zero. 

The spatial distribution of tiles in BOSS was set such that the decollided completeness of galaxy and Lya forest 
quasar samples was 100%. Due to the inhomogeneity of the target list, it is not possible to reach 100% decollided 
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completeness and 100% fiber efficiency. Since the tiling completeness was the higher priority, the BOSS fiber efficiency 
for LOWZ, CMASS, and quasar targets was ~ 90%. 

In eBOSS, we adopted a tiered-priority system for survey targets to maximize both tiling completeness and fiber 
efficiency. All non-LRG targets receive maximal priority and are assigned fibers first. The tiling solution must satisfy 
the requirement of 100% tiling completeness for the decollided set of all non-LRG targets. For LRGs, we no longer 
require full decollided completeness. Rather, the density of LRG targets intentionally oversubscribes the remaining 
fiber budget. The input sample is tiled at a lower density due to collisions with higher-priority targets, collisions 
between multiple LRGs, and lack of available fibers for LRGs in the decollided set. We refer to collisions of LRGs with 
higher-priority targets as knockouts in order to differentiate them with collisions between two LRGs. To achieve the 
survey goal of candidate LRG targets observed at 50 deg“^, the input target catalog is increased to a target density 
of 60 deg“^. Thus, after high-priority targets are assigned fibers, plates that have a surplus of unused fibers can 
sometimes be used to observe LRGs at a density higher than 50 deg“^. The quantity of interest for LRG tiling is the 
completeness of targets that are both decollided (i.e., with respect to other LRGs) and no-knockout (ie, with respect 
to high-priority targets). The former effect can be corrected in clustering measurements by proper weighting of tiled 
targets. The latter is essentially a mask overlayed on the LRG portion of the survey. In eBOSS, we require that the 
completeness of decollided, no-knockout LRGs be at least 95%. This layered tiling scheme allows 8% more area to be 
covered than would be possible if the eBOSS fiber efficiency were the same as in BOSS. 

In the first round of fiber assignments—the non-LRG targets—we must specify the priority for fiber assignments 
within collision groups. Because the quasar targets are significantly higher density than TDSS and SPIDERS targets, 
quasar-TDSS/SPIDERS collisions are fractionally more common for TDSS/SPIDERS target classes. Thus collisions 
are resolved in the following order (highest to lowest priority): SPIDERS, TDSS, reobservation of known quasars, 
clustering quasars, and variability-selected quasars. Quasars found in the FIRST survey and white dwarf stars that 
can be used as possible calibration standards are given the final priorities for resolving fiber collisions. Because the 
density of clustering quasar targets is comparable to the galaxy sample in BOSS, and because the fiber assignments 
require 100% completeness on the decollided sample, the resulting sample of clustering quasars follows tiling statistics 
that are very similar to the galaxies in BOSS. 

Five distinct regions of sky (denoted ebossN, where N is a number ranging from 1-5) were tiled in the first year 
of eBOSS. The area and tile centers for those regions are shown in Figure 3. These regions convey the average tiling 
statistics that can be expected for the quasar samples over the duration of the program. The statistics of the quasar 
samples for each region are presented in Table 3. The input target density for the QS0_CDRE sample of quasars ranges 

between 73 and 88 deg“^, a bit lower than the expected density of 90 deg“^. The density of previously-observed 

objects that satisfies the primary quasar selection is 17, 32, 31, 41, and 32 deg“^ for eboss[l-5], respectively. Region 
ebossl is outside the final SDSS-II footprint, so fewer SDSS spectra exist and the number of known targets is reduced 
by a factor of two. Generally, the variation in the density of known objects explains the variation in the new quasar 
sample; the total QS0_C0RE density only ranges from 113.1 to 118.2 deg“^ for eboss[2-5]. The total QS0_C0RE density 

in ebossl is 100.9 deg“^, significantly lower than the average of the other chunks. Regression tests indicate that this 

suppression in surface density is likely due to the higher extinction in ebossl relative to the rest of the SDSS imaging 
area. See Myers et al. (2015) for a full discussion of non-uniformity in the quasar target selection. 



Figure 3. Field centers for eboss[l—5]. The SEQUELS area is clearly defined by white space between the boundaries of eboss4 and eboss5. 
The area covered here is the area that was tiled in the beginning of SDSS-IV and the approximate survey area expected to be completed 
in the first two years of observation. 


On average, the input target sample of clustering quasars is assigned fibers at a completeness of 95%, thus reducing 
the expected number density of 0.9 < z < 2.2 quasars by 5% relative to the parent sample presented in Table 1. A 
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Table 3 

Fiber assignment statistics in the first five tiled regions for the quasar programs. Surface densities are all presented in units of deg~^. 
The “target density” corresponds to the density of objects from each class that satisfy the photometric selection criteria. The “fiber 
density” corresponds to the density of targets that were assigned fibers in the tiling process. 


Chunk 

Number of 
Plates 

Area 

(deg2) 

Total CORE^ 
Target Density 

New CORE*’ 
Target Density 

New CORE*’ 
Fiber Density 

Reobserved Lya 
Fiber Density 

PTF*^ Target 
Density 

PTF^^ Fiber 
Density 

ebossl 

199 

987.9 

100.9 

84.10 

80.47 

6.19 

18.59 

17.23 

eboss2 

98 

469.4 

113.1 

80.84 

77.08 

9.47 

10.56 

9.67 

ebossS 

204 

984.7 

118.2 

87.52 

83.69 

8.55 

21.24 

19.80 

eboss4 

80 

397.1 

114.4 

73.47 

69.86 

9.16 

28.70 

26.36 

eboss5 

70 

353.7 

116.4 

84.19 

79.68 

8.49 

14.51 

13.30 


^The sample denoted “Total CORE” corresponds to all targets satisfying the QSO.CORE selection criteria, regardless of prior observations. 

^The sample denoted “New CORE” corresponds to all targets satisfying the QSO.CORE selection criteria that were not observed previously 
and are candidates for fiber assignment in eBOSS. 

^The sample denoted “PTF” corresponds to all targets satisfying the variability selection criteria using PTF imaging data that are 
exclusive of all other samples. 

similar trend is found for the PTF-selected quasars, but at slightly lower completeness (~92%) due to their lower 
priority in resolving fiber collisions. Note that the repeated observations of Lya quasars are by definition exclusive 
of the QS0_C0RE selection but the PTF-selected quasars do have significant overlap with the QS0_C0RE sample. In 
Table 1 and Table 3, we have assigned PTF-selected targets only to the QS0_CDRE sample when overlap occurs. The 
PTF-selected densities presented in the table therefore reflect the density of unique PTF-selected targets. Because the 
FIRST sample only produces quasars at 1 deg“^, we do not include the statistics of that sample in Table 3. 

The area of the first five regions that is covered by only a single tile is 2054 deg^, 64% of the total 3193 deg^ currently 
tiled. In these areas, there is no way to capture objects lost to fiber collisions with other targets. The remaining 1139 
deg^ is covered by two or more tiles. In these regions, the completeness of the collided objects increases significantly, 
leading to a total completeness on the quasar sample of close to 100%. 

In the second round of tiling, LRGs are the only target species and require no prioritization to resolve fiber collisions. 
The statistics of the LRG sample are presented in Table 4. Because the redshift distribution of the LRG sample only 
overlaps ~ 8.5 deg“^ new QSO.CORE objects at redshifts z < 1, the masked region is mostly uncorrelated with clustering 
in the LRG sample. There is negligible overlap between the LRG and TDSS/SPIDERS samples. To first order, the 
areas of sky restricted from observing LRG targets can be treated in a similar manner to regions lost to bright stars or 
imaging artifacts in BOSS. Quantifying the full consequences of the overlapping samples will be done on mock catalogs 
as part of future clustering analyses. 

On average, the total completeness of the LRG sample is 87% while the completeness of decollided, noknock (no¬ 
knockout) targets is 98%. The average density of LRG targets given a fiber is 52 deg“^, slightly exceeding the 
goal of 50 deg“^. The resulting fiber assignments occasionally produce local fluctuations that may have insufficient 
completeness for clustering analysis. An example of the fluctuations of the decollided, noknock LRG targets tiled in 
eboss3 is shown in the left hand panel of Figure 4. The area-weighted cumulative histogram of completeness over 
the entire region covered by the first five regions is shown in the right hand panel of Figure 4. Assuming these five 
regions are representative of the completeness we expect in the full eBOSS footprint, 5% of the area will fall below 
85% completeness in the decollided, noknock LRG targets. These areas exceed the 15% uniformity requirement that 
we have generally assumed and may be excised from the clustering analysis depending on the results of future tests 
on mock catalogs. 


Table 4 

Fiber assignment statistics in the first five regions tiled for eBOSS. Surface densities are all presented in units of deg“^. 


Chunk 

LRG Target 
Density 

LRG Knockout 
Density 

LRG noknock, collided 
Density 

LRG Fiber 
Density 

Noknock, decollided 
Completeness 

ebossl 

64.66 

4.41 

3.56 

56.64 

98.9% 

eboss2 

58.73 

3.87 

2.53 

52.90 

99.7% 

ebossS 

60.53 

4.11 

2.75 

52.74 

97.3% 

eboss4 

57.12 

4.07 

2.48 

49.24 

96.6% 

eboss5 

59.37 

4.81 

2.66 

50.37 

96.4% 


In the 300-plate ELG program, all ELG targets will be assigned equal priority and ranked ahead of other targets. 
Because they have not been tiled over a cosmologically useful volume, we do not present the statistics of the ELG 
sample. Instead, those statistics will be included in a future paper dedicated to the final ELG target selection algorithm. 

4.4. Plate Design 

In the tiling stage described above, we assign science targets to fibers based on the input from target selection 
algorithms and the available fiber budget. In plate design, we convert those positions from celestial coordinates to 
the coordinate system of the telescope focal plane. We also assign additional targets to the reserved fibers to provide 
reference for sky subtraction and flux calibration. 
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Figure 4. Left: The completeness of the noknock, decollided LRG sample over the ebossS region. Right: The cumulative distribution 
of completeness in the noknock, decollided LRG sample. The distribution is weighted by the area of each independent sector defined by 
areas covered by overlapping and unique tiles. 

Tests in BOSS reveal a measurable increase in sky-subtraction residuals when the number of sky fibers drops below 
80 per plate. We therefore maintain the density of at least 80 sky fibers per plate in eBOSS. In addition, we now 
require that at least 30 sky fibers are assigned to each spectrograph. 

Standard stars for flux calibration in BOSS were selected at a density of 20 per plate over the magnitude range 
15 < rfib 2 < 19. We maintain a density of 20 F-stars selected as standard stars, but now require that at least six stars 
be assigned to both spectrographs. The F-stars at the faint limit in BOSS produced spectra that were of marginal 
quality for flux calibration, while the brightest stars had nearly three orders of magnitude higher flux than the faintest 
science targets and increase the risk of contamination from cross-talk between neighboring fibers. For this reason, 
we choose eBOSS F-stars to have a narrower magnitude distribution, with 16 < rfib 2 < 18. To compensate for the 
reduced density in candidate F-stars, we increase the scaled distance in extinction-corrected color space from the color 
of a fiducial F star by a factor of two relative to BOSS. We now require that mdist < 0.16, where 


^dist — 



5 ) - 0.82)2+ ((g-r)-0.30)' 


-b ((r -i)- 0.09)2 q_q 2)2 
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In BOSS, after accounting for predicted atmospheric differential refraction (ADR) for each plate, galaxy targets 
were centered in the focal plane to maximize throughput for 5400 A light. The hole positions for quasar targets were 
centered to maximize throughput at 4000 A to increase the S/N in the Lya forest. In eBOSS, the QS0_C0RE sample 
is selected to lie at z > 0.9 to provide direct clustering measurements while the other quasar samples are selected to 
lie at z > 2.1 to increase the Lya forest signal with respect to BOSS. All objects selected in the QS0_CDRE sample 
are centered in the focal plane at a position corresponding to the focus of 5400 A light. These targets will have 
appropriate flux calibration derived from the F-stars. All z > 2.1 quasars known from BOSS and all quasar candidates 
selected by variability are centered in the focal plane at a position corresponding to the focus of 4000 A light. While 
not appropriate for flux calibration in the current data reduction pipeline, the routines developed in Margala et al. 
(2015) can be applied to the Lya target spectra to improve the broadband distortion introduced by the ADR offset. 
As before, the wavelength that determines the center of the hole position is recorded in the quantity LAMBDA_EFF. In 
BOSS, washers manufactured with an adhesive were applied to the back of the plates at the location of quasar targets. 
Washer thickness varied according to distance from the plate center to account for the 0-300 micron difference in focus 
between 4000 A light and 5400 A light. These washers matched the quasar fibers to the focal plane for 4000 A light 
for optimal focus by offsetting the fiber tip in the direction perpendicular to the surface of the plate. However, due 
to weather, guiding, and other stochastic processes, we were unable to confirm the ability of the sticky washers to 
signihcantly improve the S/N in the Lya forest. Those washers are therefore not applied in eBOSS. 


4.5. Observing Sequence 

Simple data reductions are performed in real time to provide quick feedback to the observers and to track the 
depth of each exposure. The observers acquire signal on each field until the accumulated depth exceeds an empirically 
derived threshold for each camera. In BOSS, the depth was tuned to reach desired redshift success rates on the highest 
redshift galaxy targets in the shortest exposure time possible. Exposures were tuned so that the typical square of the 
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signal-to-noise per pixel was at least 22 for an object with native SDSS magnitude igb 2 = 21 over the wavelengths 
covered by the t-band filter. We also required the square of the signal-to-noise per pixel to be at least 10 for an object 
with (/fib 2 = 22 over the wavelengths covered by the 5 -band filter. In eBOSS, we will maintain the same thresholds 
for the beginning of the survey. To ensure that we define plate completion criteria so that we complete the goal of 
1800 plates, we will use the first year of observation to empirically determine new exposure thresholds. As was done 
in BOSS, the depths will be tuned to ensure that we can reliably classify targets in the shortest exposure possible. It 
is likely that different thresholds will be used for the quasar and LRG plates than those used for the ELG plates. 

Because quasars and LRGs are the primary focus of the eBOSS program, the majority of fibers and survey area 
are dedicated to these two target classes. The total area covered by these samples will be 7500 deg^ divided over one 
contiguous region in the SGC and one contiguous region in the NGC. The ebossl region presented in Figure 3 was tiled 
to cover an SGC area at high declination that is easily observed. The eboss2 region was chosen to obtain spectroscopy 
in a 500 deg^ region that overlaps with the Dark Energy Survey (DES) footprint. The NGC is tiled in one contiguous 
region at declination below 60 degrees. 

As was explained in Section 4.2, the sample selection for the quasar and LRG samples is complete. The first two 
years of observation will be dedicated almost exclusively to these primary fields. After two years, approximately 600 
plates should be completed for the LRG and quasar targets, producing a sample comparable in area to the DR9 
CMASS cosmology sample. Under the assumption of average weather conditions, we expect to cover almost the entire 
area covered by eboss[l-5] and shown in Figure 3. 

The third and fourth years will be split evenly between observations of ELG plates and observations of the LRG 
and quasar plates. The scheduled time in which the SGC is visible will be dedicated to ELG targets while the NGC 
time will be dedicated to the LRG and quasar targets. At a rate of 300 plates per year, we expect the ELG program 
to be completed by the summer of 2018. Because the selection of ELG targets remains uncertain, we do not define 
the exact footprint of that survey at this time. For 2016-2018 LRG and quasar observations, we expect 300 plates to 
be completed in a 1500 deg^ area of the NGC that lies just below the eboss3 and eboss5 regions. The exact area will 
depend on how much progress is made in the 2014-2015 and 2015-2016 observing seasons. 

Observations of ELG targets will be completed before the summer of 2018, regardless of whether we meet the goal 
of 300 plates. By completing the ELG program in the fourth year, we establish a meaningful cosmological sample in 
the shortest possible time. The final two years will be dedicated entirely to the LRG and quasar targets, with the goal 
of covering as much SGC and NGC area as possible. 

5. SPECTROSCOPIC DATA QUALITY 

As explained in the appendix of Alam et al. (2015a), good fortune during 2010-2014 allowed the main BOSS program 
to finish early, leaving time for new dedicated spectroscopic observations. 128 plates were drilled for the SEQUELS 
program, covering the 466 deg^ region bounded by eboss4 and eboss5 in Figure 3. In this section, we use the public 
SEQUELS data to present the quality of the LRG, clustering quasar, and Lya forest quasar spectra expected in 
eBOSS. A mosaic of randomly selected spectra spanning the relevant range of redshift for each of these three samples 
is presented in Figure 5. 

SEQUELS was designed as a pilot survey for eBOSS, using a slightly broader selection for LRG, clustering quasars, 
and Lya forest quasars that was guaranteed to include the final eBOSS selection for these classes. 66 SEQUELS plates 
were completed before the end of SDSS-III and are included in DR12. Details for accessing these data are online® and 
in Alam et al. (2015a). The remaining 62 plates will be observed at highest priority during the first year of eBOSS. 
The 466 deg^ area will not be re-tiled with eBOSS-selected targets; instead, the SEQUELS targets will be integrated 
into future clustering analysis according to the same selection algorithms as those in eBOSS. 

All SEQUELS targets are tracked by the dedicated EB0SS_TARGET0 bitmask. In the 66 plates that were observed, 
fibers were placed on 11,687 unique targets from a modified version of the final eBOSS LRG selection algorithm. 
Although the LRG selection in SEQUELS is broader than in eBOSS, the eBOSS targets can be identified by objects 
with LRG_RIW (corresponding to bit number 2), extinction-corrected magnitude imodei < 21.8, and extinction-corrected 
magnitude lUl < 17.6 (approximately 20.3 in AB). 10,873 of the SEQUELS LRG targets satisfy the final eBOSS 
LRG selection algorithm. Likewise, 19,461 unique clustering quasar targets, 6,479 PTE-selected quasar targets, and 
1,368 reobserved Lya quasars are found in SEQUELS. The clustering sample and the reobserved Lya quasar sample 
are selected in an identical fashion to those in eBOSS, identified by QS0_EB0SS_C0RE (bit 10) and QS0_RE0BS (bit 12), 
respectively. The PTE sample contains a higher density of objects than in eBOSS and is difficult to reproduce without 
the variability parameters for each object. The variability parameters will be found in the final BOSS quasar catalog 
(Paris et al. 2016) when it is complete. 

The SEQUELS spectra are the source of the redshift distributions for quasars presented in Table 1. Only SEQUELS 
data are presented in this section for the quasar samples. Both SEQUELS and proprietary eBOSS spectra were used to 
determine the redshift distributions for LRGs presented in Table 1, although the discussion of pipeline performance and 
the example spectra are based only on SEQUELS. The spectra of a large number of these objects were visually inspected 
to help characterize pipeline performance and to settle on the final redshifts of the samples. Pilot observations and 
visual inspection during the first few months of SDSS-IV led to the estimates of ELG redshift distributions presented 
in Table 2. Since the final ELG program is not decided, we postpone further discussion of those samples until future 


http: //www.sdss.org/drl2 / 
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Figure 5. Examples of SEQUELS spectra that span the range of redshifts expected in the LRG, clustering quasar, and Lya forest quasar 
samples. In each, the data are represented in black, the flux errors on each pixel in red, and the template in blue. A boxcar smoothing kernel 
of width 5 pixels has been applied to the data for illustrative effect. Each spectrum is classified with high confidence by the automated data 
reduction pipeline. Top Left: An LRG at z = 0.64. Top Right: An LRG at 2 : = 0.88. Middle Left: A quasar at 2 : = 1.08 identified 
by the QS0_C0RE selection algorithm. Middle Right: A quasar at 2 = 1.74 identified by the QS0_CQRE selection algorithm. Bottom Left: 
A quasar at 2 = 2.21 identified by the QSO.CORE selection algorithm. Bottom Right: A quasar at 2 = 3.15 identified by variability in the 
PTF imaging data. 

publications. Below, we describe the typical quality of LRG and quasar spectra and the interpretation of the automated 
redshifts and visual inspections that led to the JV(z) estimates assumed for eBOSS. 

5.1. LRG Spectra from SEQUELS 

The SEQUELS data are reduced by the same “idlspec2d” routines as those used in BOSS. The spectral templates 
described in Bolton et al. (2012) are fit to each reduced spectrum to derive a redshift and object classification. These 
templates were tuned for optimal performance in the BOSS galaxy sample and exceeded the predicted redshift success 
rates for even the most distant galaxies. The automated classification of the final BOSS CMASS galaxy sample shows 
98.4% of objects with ZWARNING_N0QS0 equal to zero, indicating that the automated redshift estimate is reliable when 
using a template suite that excludes quasar templates. 

The choice of spectral templates was sufficient for BOSS, but is not optimized for the fainter, higher redshift LRG 
galaxies that comprise the eBOSS LRG sample. Of the 10,873 eBOSS LRG targets in SEQUELS, 7,796 produced 
a ZWARNING_NDQS0 value equal to zero in the automated reductions. A statistical error on the redshift estimate is 
provided by the automated classification routine for each of these galaxies. Tests in BOSS using repeat observations 
showed that these errors were underestimated by up to a factor of 1.34 for galaxies (Bolton et al. 2012). Assuming 
that the SEQUELS/eBOSS errors are underestimated by a factor of 1.34, we find that the median redshift error is 
only 58 km s“^ and that only ten objects have a redshift error larger than 200 km s“^. The automated classification 
characterizes redshifts with a precision much better than required. 

The single largest failure mode for the automated classifications is due to the inability to discriminate between best 
fit templates at different redshifts. Objects that fail according to this criteria are assigned a warning flag equal to two, 
which means the difference between the best spectral fit and the second best spectral fit is less than 0.01 per degree 
of freedom. This failure mode is not unexpected; the eBOSS LRG targets are intrinsically fainter than the BOSS 
galaxies and thus have a lower S/N measurement of flux at each pixel in the spectrum. The faint spectra result in less 
pronounced absorption features that would otherwise provide a good redshift estimate and a reliance on broad band 
flux that is harder to fit as a function of redshift. The overall reduction of the difference with respect to the BOSS 
galaxy sample is shown in the left panel of Figure 6. The full distribution clearly shifts toward lower Ay^, leading to 
a higher fraction of objects failing to meet the required difference of 0.01 per degree of freedom. 

Of the 7,796 objects that were given a ZWARNINGJIDQSO value equal to zero, only 6,096 fall in the redshift range 
0.6 < z < 1.0. The corresponding 56.1% combined targeting and spectroscopic efficiency would fall below our 
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Figure 6. Left: The distribution of Ax^ per degree of freedom for the BOSS sample of CMASS galaxies and the eBOSS LRG sample. 
The x^ per degree of freedom value represents the difference between the best fit spectral template and the second best spectral template 
when quasar templates are excluded from the fits. Right: A scatter plot showing the difference in redshift {kms~^) between pairs of 
observations taken of the same LRG target. The two vertical lines represent threshold for the current Ax^ ZWARNINGJIOQSO (blue) and the 
more lenient threshold that increases the rate of redshifts that are considered reliable (red). The horizontal line (magenta) represents the 
limit at which redshift discrepancy is considered a catastrophic failure. 

requirement (80%) to obtain 300,000 spectroscopically-confirmed LRG tracers. In particular, the fact that 1700 of the 
10,873 spectra are given ZWARNING_NQQS0 equal to zero and do not lie in the redshift range of interest makes it nearly 
impossible to meet the requirement for 80% completeness. 

We further explored the quality of the automated classification in two ways to see how close the selection is to the 
required tracer density. First, we identified 591 LRG targets that were tiled on more than one plate and produced 
multiple spectra. We supplemented this sample by re-running the data reduction pipeline on four plates that were 
observed to more than twice the normal depth, dividing those exposures into unique coadds of the data. The data split 
produced 630 targets with multiple spectral classifications. We assessed the rate of catastrophic redshift failures by 
counting the fraction of objects for which a pair of observations produced discrepant redshifts. For the sample of 2442 
total observations, we found 166 discrepant redshifts, corresponding to a catastrophic failure rate of 6 . 8 % if we assume 
that consistent redshifts are correct. We further explored the rate of discrepant redshifts as a function of per degree 
of freedom. As shown in the right panel of Figure 6 , the vast majority of discrepant redshifts occur for rchisq_N0QSD 
< 0.005, implying that the threshold of 0.01 may be too strict. Gonsidering that a spectrum has 4400 pixels over 
the range 3600 < A < 10000 A before accounting for masked artifacts, reducing the threshold to 0.005 corresponds 
roughly to = 20. Filtering on targets with this new threshold, we find that the rate of ZWARNING_N0QSD failures 
drops from 28% to 16%. Of the 1650 observations of objects where both repeated spectra were fit with rchisq_N0QSD 
> 0.005, only ten produced discrepant redshifts, implying a catastrophic failure rate of 0.6%, below the maximum 
allowed ( 1 %). 

When reducing the threshold for rchisq_N0QS0, we reduce the rate of ZWARNING_N0QSD failures by nearly a factor 
of two and appear to meet the requirements for catastrophic redshift failures. However, even with this change we are 
still well below the required rate of spectroscopic confirmation. As a second test, we proceeded to visually inspect the 
LRG spectra from 18 different plates that obtained the deepest observations. The visual inspections were intended to 
provide insight into the spectral quality, the properties of the templates producing the best and second best fits, and 
to gauge the true redshifts of the LRG sample. 

While it is impossible to run visual inspections in a reproducible and consistent fashion, we attempted to provide a 
scheme by which we could assess the range of likelihoods that a redshift is correct. We established a simple four-tiered 
redshift confidence metric Zconf- 

• -s^conf = 0 denotes a spectrum for which it is impossible to classify the type or determine a redshift. Of the 
LRG target spectra that were visually inspected, 6 . 6 % were assigned Zconf = 0. 88 % of these spectra were given 
ZWARNING_N0QS0 > 0 in the automated classification. 

• -Zconf = 1 denotes a spectrum that is ambiguously classified, where more than one of the best-fit models is a 
possible fit. The redshift provided by the inspector is intended to be a best estimate of the correct redshift. 
4.5% of the spectra were determined to be galaxies and assigned Zconf = 1- 48% of these spectra were given 
ZWARNING_N0QSD > 0 in the automated classification. 

• -Zconf = 2 denotes a spectrum where the redshift estimate is likely to be correct as there is no other template that 
provides an equally compelling fit. However, some degree of uncertainty remains because the best-fit spectral 
template does not reproduce all of the expected features. For example, the template may fit the continuum but 
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certain absorption features may be ambiguous in the data. 12% of the spectra were determined to be galaxies 
and assigned Zconf = 2. 23% of these spectra were given ZWARNING_NOQSD > 0 in the automated classification. 

• Finally, Zconf = 3 denotes a case where the redshift is estimated at very high confidence. 68% of the spectra were 
determined to be galaxies and assigned ^conf = 3. 4.8% of these spectra were given ZWARNING_N0QS0 > 0 in the 
automated classification. 

• 8.9% of the objects were classified as stars. M-stars were the dominant stellar contaminant and were very easy 
to identify in visual inspection. 92% of the stars were given Zconf = 3, but 42% were given ZWARNING_NOQSD > 0. 
The high rate of ZWARNING_NOQSD was later found to be caused by a bug in the version of the pipeline used for 
DR12. This will be fixed for the next public data release. 


Two redshift distributions are presented in Table 1. The more conservative estimate (higher rate of “Poor Spectra”) 
assumes that any target given Zconf > 1 is given the correct redshift. The less conservative estimate corresponds to all 
objects with Zconf > 0. It is likely that the true distribution lies between these two bounds. 



4000 5000 6000 7000 8000 

Observed Wavelength (A) 



4000 5000 6000 7000 8000 

Observed Wovelength (A) 


Figure 7. Example of the influence on classification by non-physical flux at extrema of spectrum. The data are represented in black, the 
flux errors on each pixel in red, and the template in blue. A boxcar smoothing kernel of width 11 pixels has been applied to the data for 
illustrative effect. Left: The template corresponding to the visual inspection redshift is fit to the data with a degree of freedom of 

1.017 at a redshift z = 0.80831. Right: An incorrect template is fit to the data with degree of freedom of 1.022 at ^ = 1.09220, a 

difference from the first template small enough to trigger the ZWARNING_N0QS0 flag. 


Visual inspection reveals at least three failure modes among the targets with small Ay^. One failure mode occurs 
when the spectrum possesses a non-physical gradient in the flux at the bluest wavelengths or the reddest wavelengths. 
An example of one such failure is shown in Figure 7. One notes a clear identification of the Ca H&K features around 
observed wavelength 7200 A in the left hand panel. However, the data deviate significantly from the template over 
the range 3600-4000 A, leading to an inflated estimate. The template in the right hand panel does not capture 
the narrow absorption features because the redshift is incorrect. Instead, the template is better fit to the excessive 
UV flux that is likely due to imperfect spectral extraction and not of astrophysical origins. In cases such as this, the 
polynomial term that accounts for uncertainties in flux calibration appears to be better coupled to incorrect galaxy 
templates than to a galaxy template at the correct redshift. 

A second failure mode is shown in Figure 8. Because the sources in the LRG target sample are typically passive 
galaxies at a high redshift, the significance of the 4000 A break can be diminished due to low flux counts. In these 
cases, it becomes difficult to discriminate between the 4000 A break and G-band absorption extending to 4304 A. A 
spectral template of a higher redshift galaxy with a suppressed G-band feature (right panel of Figure 8) can sometimes 
fit the observed spectrum nearly as well as a template with a stronger G-band feature and a correct fit to the 4000 
A break. While the template that produces the stronger G-band absorption profile is more physically motivated than 
the higher redshift template, the current reduction pipeline has no mechanism in place to impose a prior in favor of 
the correct redshift. 

Finally, the third failure mode occurs because the majority of the signal is found in the broad band flux as opposed 
to the monochromatic Ga H&K features. In cases such as the spectrum shown in Figure 9, there is little power to 
statistically discriminate between a fit that correctly places both Ga H&K features and a fit that staggers the template 
with an incorrect fit to a single feature. In many cases, the G-band absorption line (4304 A), Magnesium line (5175 
A), and Sodium line (5894 A) provide additional constraints on the redshift, but these lines become ambiguous at 
higher redshifts where sky subtraction residuals are more significant. 

As will be discussed in Section 6.2, we are now improving the algorithms for spectral extraction and classification to 
mitigate these classes of failure modes. We are confident that we can increase the rate of reliable classifications beyond 
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Figure 8. Example of confusion between broad G-band absorption and the 4000 A break. The data are represented in black, the flux 
errors on each pixel in red, and the template in blue. A boxcar smoothing kernel of width 11 pixels has been applied to the data for 
illustrative effect. Left: The template corresponding to the visual inspection redshift is fit to the data with a per degree of freedom of 
1.258 at a redshift z = 0.64375. Right: An incorrect template is fit to the data with P®r degree of freedom of 1.264 at 2 : = 0.82202, a 
difference from the first template small enough to trigger the ZWARNING_N0QS0 flag. 
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Figure 9. Example of the spectra in which the templates are unable to differentiate the Ca H&K features at high significance. The data 
are represented in black, the fiux errors on each pixel in red, and the template in blue. A boxcar smoothing kernel of width 11 pixels has 
been applied to the data for illustrative effect. Left: The template corresponding to the visual inspection redshift is fit to the data with a 
X^ per degree of freedom of 0.911 at a redshift z = 0.61516. Right: An incorrect template is fit to the data with x^ per degree of freedom 
of 0.915 at 2 : = 0.62806, a difference from the first template small enough to trigger the ZWARNING_N0QSD flag. 

what is presented in DR12 and are optimistic that we can achieve statistics comparable to the Zconf > 0 redshift 
identification rate from the visual inspections. 

5 . 2 . Quasar Spectra from SEQUELS 

Quasar targets in BOSS were observed at a rate of roughly 200 per plate, totaling nearly 500,000 over the full 
program. Each of these spectra was visually inspected and given a classification (quasar, star, or galaxy) and a 
redshift estimate with a documented degree of confidence (Paris et al. 2012, 2014). We also relied on visual inspection 
to flag quasar sightlines with peculiarities such as damped lyman-alpha and BAL systems that could not be classified 
with the automated redshift classification (Bolton et al. 2012). The photometric and spectroscopic parameters of each 
object were recorded for use by all members of the collaboration and released in public form on a regular basis. The 
process was undertaken primarily by two members of the collaboration and proved to be very time-consuming. 

The quasar redshift distributions presented in Section 4.2 were derived from visual inspections of all SEQUELS 
plates observed during BOSS. Even the faintest quasars (rpsp = 22) can be confidently classified in visual inspections; 
95.7% of the full quasar sample (including TDSS, SPIDERS, and PTE-selected objects) can be identified with high 
conhdence. A comparison of the automated pipeline to the results of visual inspection reveal a high level of agreement. 
The redshift estimates from the spectroscopic pipeline are consistent with the visual inspection estimates to within 
1000 kms“^ for 98.7% of QS0_C0RE quasars (rpsF < 22). At this early stage, the spectroscopic pipeline nearly meets 
the requirement presented in Section 3.2 that < 1% of objects are given catastrophically incorrect redshifts. 

The density in eBOSS and SEQUELS increases to roughly 600 quasar targets per plate which motivated the collab¬ 
oration to revisit the automated classification schemes and significantly decrease the level of visual inspection. To fully 
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transition to this new scheme, we evaluated the pipeline performance against visual inspections to determine the source 
of catastrophic failures and systematic uncertainties in redshift estimates. We found patterns in catastrophic failures 
associated with the class of object reported by the pipeline, the redshift estimate, and the second, third, fourth and 
fifth best estimates. By applying the following criteria, we are able to improve the consistency between the pipeline 
classifications and the visual inspections. 

• Objects identified as a star in the best estimate are assumed to be a star. 

• Objects identified as a galaxy with 2 ; < 1 in the best estimate are assumed to be a galaxy at 2 ; < 1. 

• Objects identified as a galaxy with 2 ; > 1 in the best estimate are assumed to be a galaxy if one of the next four 

classifications is also a galaxy. 

• Objects identified as a quasar with ZWARNING= 0 are assumed to be a quasar unless two or more alternate 
identifications are given a class of star. In these cases, the object is assumed to be a star. 

• Objects identified as a quasar with ZWARNING> 0 are assumed to be a star if two or more alternate identifi¬ 
cations are given a class of star. 

• Objects that fail to meet these criteria are followed up with visual inspection. 

The results of this classification scheme applied to SEQUELS data is presented in Table 5. This new classification 
provides identifications for 93% of all the quasar targets. Using the visual inspections as a baseline, among those auto¬ 
matically classified targets, 96.9% are correctly characterized by our new algorithm. It is specifically designed to pick 
up actual quasars and to limit the fraction of contaminants in the quasar sample and lost quasars: 54.7% are actual 
quasars that are confirmed after visual inspection. A total of 0.42% quasar targets are wrongly classified as quasars 
by our automated scheme. Those contaminants are largely dominated by spectra with very low signal-to-noise ratios, 
and could also not be identified after visual inspection. Our procedure misidentifies 0.40% of quasar targets that are 
confirmed to be quasars after visual inspection. Among the 7% of all quasar targets that still require visual inspection, 
70.3% are true quasars. Assuming that the identification error from the visual inspection is negligible, the total 
fraction of contamination of the quasar sample and loss is below the 1% of the overall sample as required in Section 3.2. 


Table 5 

Performance of the automated classification scheme applied to the SEQUELS quasar sample. Percentages are computed with respect to 
the total number of automatically identified objects for the columns STAR, QSO and GALAXY. For objects without an automated 
classification (UNKNOWN), percentages are calculated with respect to the total number of UNKNOWN objects. Quasar targets 
classified as UNKNOWN by the automated classification are visually inspected. 


Automatic classification 



STAR 

QSO 

GALAXY 

UNKNOWN 

Visual inspection result 
?/BAD 

78 (0.2%) 

102 (0.3%) 

427 (1.3%) 

486 (19.9%) 

STAR 

9,845 (29.2%) 

19 (0.05%) 

159 (0.5%) 

189 (7.7%) 

QSO 

37 (0.1%) 

18,475 (54.7%) 

103 (0.3%) 

1,718 (70.3%) 

GALAXY 

13 (0.04%) 

24 (0.07%) 

4,390 (13.0%) 

50 (2.1%) 


So we can expect that eBOSS quasar redshifts will be robust, but an additional requirement noted in Section 3.2 
is that they should also have a precision of order 300-600 kms“^. Based on repeat spectroscopy, the RMS scatter of 
BOSS pipeline redshift estimates for quasars is a factor of 2 higher than the reported statistical errors (Bolton et al. 

2012) . Taking into account these underestimated errors, the typical statistical precision is between ay ~ 100 and 
ay ~ 300 km s“^ in the redshift range 1 < 2 : < 2 , depending on the luminosity of the source. If truly representative 
of the redshift errors, this level of precision would be sufficient to avoid degradation in the clustering signal. 

The reported statistical precision does not necessarily reflect the true accuracy on a quasar redshift due to possible 
systematic uncertainties in the redshifting templates (e.g. Hewett & Wild 2010). Highly ionized gas in the broad¬ 
line region of quasars is subject to radiation-driven winds (e.g. Proga 2003). It is therefore likely that the measured 
redshifts largely determined by these emission lines are offset from the systemic redshift. Such an effect has been 
reported in the redshift estimates of 2 : > 2.1 quasars used in Lya forest - quasar cross correlations (Font-Ribera et al. 

2013) . Because Mgll is ionized at lower energies and is largely present at larger separations from the central black 
hole, Mgll emission lines are likely the least affected broad emission lines by systematic offsets. Various studies have 
estimated that the Mgll emission line is offset by an amount varying from -100 to -1-50 kms“^ with respect to the 
quasar systemic redshift (e.g. McIntosh et al. 1999; Hewett & Wild 2010; Font-Ribera et al. 2013). The upper bound 
on the statistical errors of quasar redshift estimates can be therefore estimated directly from fits to the Mgll emission 
lines. 

In the redshift range of interest for eBOSS quasar clustering, the Mgll emission fine is always covered in the 
spectra. Hence, it can be directly used to estimate quasar redshifts. In order to estimate the statistical error on 
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Mgll-based redshifts, we use the spectra of 472 quasars in the redshift range 0.9 < z < 2.2 taken as part of the SDSS 
Reverberation Mapping Project (Shen et al. 2015). 30 spectroscopic epochs were taken in SDSS-III for each of those 
quasars with a homogeneous S/N distribution corresponding to twice the normal BOSS and eBOSS exposure depth. 
The redshift and g-band magnitude distributions are comparable with the eBOSS sample. By fitting the Mgll line 
with a linear combination of five principal components, its redshift is measured at the position of the maximum flux 
of the emission line as described in Paris et al. (2012). We estimate the RMS scatter of the Mgll-based redshifts as 
a function of g-band magnitude and quasar redshift in Figure 10 (red diamonds). The redshift errors increase with 
g-band magnitude and vary from ^80 to ^300 kms“^. This behavior is expected since the S/N decreases for fainter 
objects. Redshift errors also increase from ~130 kms“^ at z = 1 to ^270 kms“^ at z = 2.2. Larger redshift errors for 
quasars at z > 2 are measured because the Mgll line lies in the red part of eBOSS spectra where the sky subtraction 
is noisy and/or imperfect and makes the line fitting more difficult. Redshift errors are also larger at z ^ 1.2 when the 
Mgll emission line moves from the blue arm to the red arm of the eBOSS spectrograph. 
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Figure 10. The statistical error on redshift estimates for the quasar sample for template-based (black squares) and Mgll-based (red 
diamonds) redshift estimates. Left: The redshift error as a function of g-band magnitude. Right: The redshift error as a function of 
redshift. 


The statistical errors on redshift measurements based on Mgll emission line fits are well within the requirements 
described in Section 3.2. However, this emission line is not always detected, especially in the spectrum of faint quasars 
at z > 2 when the Mgll line is at its red end. We thus consider redshifts estimated with quasar templates in order to 
overcome this issue. These templates have been calibrated with the Mgll emission line (see details in Paris et al. 2012). 
We estimate the statistical error on those redshift measurements in the same manner as for the Mgll-based estimates. 
As expected, template-based redshifts are more stable (black squares in Figure 10). The RMS scatter increases from 
^60 to ^100 kms“^ for 18 < 5 < 22 (left panel). There is no obvious redshift dependence of the statistical error with 
an average of 90 kms“^ (right panel). Despite this apparent better behavior of statistical errors for template-based 
redshifts, one significant drawback comes from astrophysical effects. Intrinsic emission line shifts depend on quasar 
luminosity (e.g. Hewett & Wild 2010). Offsets can be accurately reproduced with templates if the whole luminosity 
range of interest is covered by the training sample, which is not the case of our current templates. The most affected 
emission line is CIVA1550 that is blueshifted by several hundreds of kms“^ with respect to the systemic redshift (e.g. 
Richards et al. 2011). When this line enters in the spectrum, it may affect the accuracy of template-based redshift 
measurements and thus, it may cause an additional source of systematics that is redshift dependent. In order to 
quantify this effect, we measure Mgll- and template-based redshifts of quasars in the SEQUELS sample. We then 
measure the redshift evolution of the scatter between Mgll-based and template-based redshifts (Figure 11). Between 
z = 0.9 and z = 1.5, redshift errors are below 300 kms“^. At larger redshift, the CIV line is covered by eBOSS spectra 
and redshift errors increase up to 600 kms“^. This demonstrates that template-based redshifts are dominated by 
astrophysical effects. Templates need to be improved to reduce these systematics. Nevertheless, redshift requirements 
(and high level BAO projections) noted in Section 3.2 are set to our current redshift accuracy (red line in Figure 11). 




eBOSS 


27 



Figure 11. Redshift evolution of the statistical error on systemic redshift estimates for the quasar sample based on the comparison of 
template-based and Mgll-based redshit estimates. Black points are the statistical errors derived from the data. The red line shows the 
redshift evolution used for the cosmological forecasts described in Section 3.2. 

5.3. Lya Forest Spectra from SEQUELS 

The Lya forest sample differs from the LRG, ELG, and clustering quasar samples in that it will be used to increase 
the density of a BOSS sample that is already being used to measure the BAO distance scale. The analysis tools and 
spectral classification are well advanced for this sample of targets and are described in recent results. For example, 
three techniques for determining the underlying quasar continuum in the Lya forest region were compared in the latest 
BAO measurement (Delubac et al. 2015). Blomqvist et al. (2015) explore a model in Fourier space to correct broadband 
distortion introduced by continuum fitting. Implementing the technique on simulated spectra, they find that the linear 
bias parameter bp and the redshift-space distortion parameter f3p can be recovered with systematic errors less than 
0.5%. Astrophysical effects (quasar continuum diversity and high-density absorbers) and instrumental effects (noise, 
spectral resolution, and data reduction artifacts) that affect estimates of the Lya forest auto-correlation function were 
quantified through catalogs of mock spectra (Bautista et al. 2015). Analysis in SDSS-IV will be performed jointly 
on the BOSS and eBOSS Lya forest samples; modeling of the astrophysical and instrumental contributions to the 
measured clustering signal will continue to evolve. 

The sample of z > 2.1 quasars also provided the first opportunity in SDSS to measure BAO through cross-correlation 
of different tracers. The target sample produces both direct tracers of the underlying density field (quasars themselves) 
and illumination of neutral hydrogen through the Lya forest. The large-scale cross-correlation of quasars with the Lya 
forest absorption was first measured to comoving separations of 80 /i“^Mpc (Font-Ribera et al. 2013). The technique 
was later scaled to larger separations with a larger sample of quasars to measure BAO and provide constraints on 
H{z) and ^^(z) as discussed in Section 2.2. 

The measurement of BAO with Lya forest in auto-correlation and in cross-correlation with quasars are both shot- 
noise limited and produce nearly independent constraints on the distance scale. The final BOSS DR12 sample is 
expected to produce combined constraints on H{z) and dA{z) with a precision of 2.0% and 2.5%, respectively. One 
can project the relative improvement offered by the new eBOSS z > 2.1 quasar sample by evaluating the number 
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density and depth of spectra compared to those in BOSS. The observed number density of quasars can be computed 
from the redshift distributions found in Table 1 and the tiling efficiencies in Table 3. The density of quasars from each 
sample after accounting for tiling efficiency is shown in Table 6. SEQUELS observations are used to assess the depth 
of spectra. 

We follow the formalism presented in McQuinn & White (2011) to estimate the fractional improvement on the BOSS 
BAO distance precision from the auto-correlation when adding eBOSS Lya forest spectra. The BAO uncertainties are 
forecast based on the number of modes available given the survey volume and the signal-to-noise per mode (S/N). In 
the case of Lya forest spectra, the S/N scales as 


(S/N)^ 


pm 

mtik) 


( 2 ) 


where Ppik) is the flux power spectrum. Ptot{k) is the variance of a given mode approximated by the combination 


Ptotik) = Ppik) + , 

^eff 


( 3 ) 


where Pid(A:||) = Pm is the ID power along the line of sight. The effective angular density of lines of sight, UeS, is 
defined 

Ues = < V >= n < Pid/{Pid + Pn) > ■ ( 4 ) 

The average surface density of quasars used to measure absorption in the Lya forest is represented by and Pn is 
a weighted average of the noise power {Pn = dx/{SNR)'^). The mean pixel width within the forest is recorded as dx 
(units of Mpc/h) and SNR is the mean signal-to-noise ratio per pixel within the forest. In a survey that is dominated 
by shot noise, the S/N per Fourier mode should be roughly proportional to rieff at a given redshift. Although the 
calculations are redshift dependent, we assume that the redshift distributions in eBOSS are similar to BOSS so the 
uncertainties of eBOSS relative to BOSS should scale as the inverse ratio of rieff in the areas of overlap. 

In studying the precision of BAO measurements in the Lya auto-correlation with the addition of eBOSS, we simply 
determine in the area {AgBOSs) covered by BOSS and eBOSS and in the area (Aboss) covered 

uniquely by BOSS. Since the two regions will be independent, we can compute the uncertainty on the BAO in eBOSS 
relative to BOSS according to 


2 

^eBOSS 

2 

'^BOSS 


_ Aboss {nfm? _ 

AeBOSS + ^boss - A^oss) 


1 

Ia + (1 - Ia) ’ 


( 5 ) 


where we have defined Ja = AeBOSS/A boss and /„ = 


Table 6 

Densities of quasars and lines of sight to Lya forest in BOSS and the various eBOSS samples. The density of quasars ( tiq ) includes BAL 
quasars at redshifts 2 : > 2.0. The density of lines of sight {np) include only quasars with 2 > 2.1 that are not BAL quasars. BOSS 
numbers include the true non-BAL population while the eBOSS numbers include a 13% reduction in the observed number density to 

account for the predicted rate of BAL contamination. 


Description 

Symbol 

Value 

Total area of BOSS survey 

Aboss 

10 000 deg2 

BOSS area covered after 2 years of eBOSS 

AeBOSSl 

3 000 deg^ 

BOSS area covered after 4 years of eBOSS 

AeBOSS2 

4 500 deg2 

BOSS area covered after 6 years of eBOSS 

AeBOSSS 

7 500 deg2 

Density of quasars in BOSS 

Q 

'^BOSS 

71 ^ 

'‘‘C^RE 

^ptf 

19.7 deg“2 

Density of new quasars in eBOSS from CORE selection 

10.1 deg“^ 

Density of new quasars in eBOSS from PTF selection 

2.8 deg“^ 

Density of Lya lines of sight in BOSS 

'^BOSS 

16.7 deg-^ 

Density of Lya lines of sight in BOSS not re-observed in eBOSS 

'^bright 

8.7 deg-2 

Density of Lya lines of sight in BOSS re-observed in eBOSS 

^ faint 

8.3 deg-2 

Density of new Lya lines of sight in eBOSS from CORE selection 
Density of new Lya lines of sight in eBOSS from PTF selection 

'^CpRE 

Rp'-pp 

5.5 deg“^ 

2.3 deg-2 

Mean value of a Lya quasar in BOSS 

Wboss) 

0.60 

Mean value of a Lya quasar in BOSS not re-observed in eBOSS 

{^bright) 

0.85 

Mean value of a Lya quasar in BOSS before re-observed in eBOSS 

{^before) 

0.35 

Mean value of a Lya quasar after re-observed in eBOSS 

(y after) 

0.55 

Mean value of a new Lya quasar in eBOSS from CORE selection 

{ynew ) 

0.39 

Mean value of a new Lya quasar in eBOSS from PTF selection 

{ynew ) 

0.50 

Effective density of lines of sight in BOSS 

^efT 

10.0 

Effective density of lines of sight in eBOSS (H-BOSS) 

^eBOSS 
^cff_ 

15.3 
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The depth of these SEQUELS measurements for each target sample was used to compute the average effective value 
of each Lya forest line of sight according to Equation 4. Using the observed densities found in Table 6, we find 
= 10.0, = 15.3, and /„ = 1.53. The area covered in eBOSS will progress according to the survey 

strategy presented in Section 4.5. Ja will be equal to 0.3, 0.45, and 0.75 after the first, second, and third public data 
releases, respectively. 

In the case of the cross-correlation, the S/N per Fourier mode should be roughly proportional to the sampling of the 
density field. Since both the number density of quasars and the measurements of the Lya forest contribnte, the term 
rieff from the auto-correlation gets replaced by what is effectively the geometric mean between the effective number 
density of Lya forest lines of sight and the quasar density. The last form of Equation 5 can be simply rewritten as 

^eBOSS _ Cg'l 

^BOSS fA fn /« + (!- /a) 

where we have defined the additional parameter as the density of quasars in eBOSS relative to BOSS, fq = 
^‘eBOSs/^'boss- values of /a for each incremental data release are of course the same as in the auto-correlation. 
Using the number density of quasars in the second block of Table 6, we find n^Qgg = 19.7 deg“^, n^soss ~ deg“^, 
and fq = 1.65. The additional quasars from eBOSS will therefore have a slightly higher impact on the cross-correlation 
measurement than the auto-correlation measurement. 

For the auto-correlation, we expect the error on the BAO distance scale from BOSS to be reduced by factors 1.18, 
1.27, and 1.42 for the two year, four year, and six year increments in the program. For the cross-correlation, we expect 
the error to be reduced by factors 1.21, 1.30, and 1.46. Because the improvements in the auto- and cross-correlation 
are so similar, we average the two and assume improvements of 1.19, 1.28, and 1.44 on the distance scale determined 
from the joint analysis. These values are used to project the BAO distance precision from Lya forest data in Section 7. 

6. SURVEY METRICS AND FUTURE DEVELOPMENT EFFORTS 
6 .1. Data Quality Compared to Requirements 

The high-level goals of eBOSS are to obtain BAO distance measurements with 1% precision using an LRG sample, 
with 2% precision using an ELG sample, and with 2% precision using a quasar sample. A number of low-level metrics 
must be satisfied to meet these goals. The redshift range and number of spectroscopically confirmed objects required 
for each class of tracer was presented in Section 3.1. The requirements on spectroscopic classification and uniformity 
in targeting were presented in Section 3.2. 

Significant testing of the LRG and quasar samples has allowed us to determine whether we meet all of the goals. 
A summary of the target population was outlined in Section 4.2 while the expected rate of fiber efficiencies was 
presented in Section 4.3. Combining the number density of the parent target population with the tiling statistics, we 
can predict the number of spectroscopically confirmed targets that can be used as tracers. Pilot studies described 
in Section 5 produced a sample of more than 10,000 LRG and more than 10,000 quasar spectra from which we can 
estimate the performance of spectroscopic classification. Assuming a joint LRG and quasar survey covering 7500 deg^, 
we summarize the expected performance of these two target classes with respect to requirements in Table 7. As can be 
seen, we satisfy the pipeline requirements on both the LRG and quasar samples and we exceed the required number 
of clustering quasars by 15%. We expect to miss the requirement for the number of new spectroscopically confirmed 
LRGs by 12%, leading to an expected degradation in BAO signal of roughly 6% relative to the requirement. The 
overlap with the high redshift tail of the BOSS GMASS sample increases the total number of 0.6 < z < 1.0 passive 
galaxies by 66%, thus allowing high precision measurements of BAO over this redshift range. 

The final ELG selection remains undecided, so it is not yet possible to present the expected performance for this 
target class. We do have estimates for the number density using several techniques presented in Section 4.2. We 
present the statistics we can estimate for two of those samples in Table 7. We report the statistics for the Fisher 
UgrizWl covering 1500 deg^ and the high density DECam selection covering 1100 deg^. Because the emission lines 
are narrow, we expect to easily meet the requirement of < 300kms“^ redshift precision, although we do not have a 
specific estimate at this time. More challenging will be robust identification of sources; catastrophic failures due to line 
confusion and sky subtraction residuals pose a risk to this class of target. We postpone discussion of the catastrophic 
failure rate and the expected uniformity until a future publication. 

6 .2. Future Improvements in Automated Data Reduetion 

The eBOSS spectroscopic pipeline will be based on the BOSS spectroscopic pipeline that was adapted from the orig¬ 
inal SDSS “idlspec2d” pipeline. It consists of two main parts: the “spectro2d” portion that extracts two-dimensional 
raw data into one-dimensional flux calibrated spectra; and the “spectrold” portion that measures scalar quantities 
such as classifications and redshifts from those one-dimensional spectra (Bolton et al. 2012). For the most part, eBOSS 
requires the same spectro2d performance as achieved in BOSS, but extended to fainter targets. As discussed in Sec¬ 
tion 5.1, a number of problems in extraction and/or flux calibration have arisen that degrade the redshift classification. 

New development on the data reduction pipeline has begun to address these problems. This renewed effort includes 
removing a S/N dependent bias in the extractions, coadditions, and sky subtraction; improving the sky model in the 
presence of variations in brightness and fiber resolution over the focal plane; improving the error propagation in the 
coadditions; and correcting occasional failures of flux calibration at the extreme red or blue ends of the spectrum. New 
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Table 7 

Comparison of expected eBOSS performance compared to requirements for high-level BAO measurements. 



LRG 

CMASS 

Clustering Quasars 

UgrizWl ELG 

DEGam ELG 

Redshift Range 

0.6 < 2 < 1.0 

0.6 < 2 < 1.0 

0.9 < 2 < 2.2 

0.6 < 2 < 1.0 

0.7 < 2 < 1.1 

Target Density and Fiber Density 

Density of Targets (deg“^) 

60 

23 

115 

180 

240 

Density of Previously Observed Objects 

0 

23 

25 

0 

0 

Density of Objects Assigned Fibers 

52 

0 

85 

171 

228 

Efficiency of Redshift Classification 

Density Expected Confirmed Tracers (deg~^) 

36 

23 

67 

121 

179 

Number of Required Confirmed Tracers 

300,000 

- 

435,000 

190,000 

190,000 

Number of Expected Confirmed Tracers 

266,000 

174,000 

500,000 

182,000 

197,000 

Statistical Accuracy of Redshift Estimates 

Required Redshift Accuracy (kms 

< 300 

< 300 

< 300-I-400[2- 1.5] 

< 300 

< 300 

Expected Redshift Accuracy (kms~^) 

58 

35 

< 300 4-400(2- 1.5] 

< 300 

< 300 

Unrecognized Redshift Classification Failures 

Catastrophic Failures Requirement 

< 1% 

< 1% 

< 1% 

< 1% 

< 1% 

Catastrophic Failures Expected 

0.6% 

< 1% 

0.8% 

TBD 

TBD 

Uniformity Over Targeting Area 

Uniform Area in Imaging Systematics 

92% 

- 

90% 

TBD 

TBD 

An/A(0.01ma^) Zeropt error 

6.2% 

- 

0.86% 

TBD 

TBD 

Uniform Area in zeropoint errors 

86.7% 

- 

100% 

TBD 

TBD 


flux calibration that accounts for ADR across the focal plane (Margala et al. 2015) will also be implemented in the 
data reduction pipeline. These changes will provide more accurate spectra with better error models. Without making 
any changes to the spectral templates, we hope to resolve a significant fraction of uncertain spectral classifications for 
the new faint targets in eBOSS by making these changes to the spectro2d pipeline. 

The analysis of quasar Lyman-alpha forest beyond the baseline BAO measurements may require further developments 
to the spectro2d pipeline, though the exact requirements are not specified at this point in time. Characterization of 
the spectra and development of new algorithms falls under the purview of the Lyman-alpha working group and not 
the team that develops the data reduction pipeline. 

Any new developments in the spectro2d component of the data reduction pipeline will be backward compatible with 
BOSS data, enabling consistent reprocessing for joint analyses of BOSS and eBOSS data. We do not require SDSS-I 
and SDSS-II data to be re-processed because the smaller wavelength range renders the early z > 2.15 quasar spectra 
obsolete with respect to the BOSS and eBOSS spectra. Because of a lack of refined selection techniques and reduced 
system throughput, no high redshift ELG or LRG spectra are found in the earlier programs. 

The spectroscopic classifications and redshifts from spectrold provide the final map for the clustering measurements 
and must meet the requirements laid out in Section 3.2. As described throughout the text and in Table 7, the BOSS 
spectral templates are nearly sufficient to meet these goals. BOSS quasar identifications were visually checked by 
two individuals to unambiguously classify the BOSS quasars to a level of detail not currently possible by automated 
software. There is a small sample of QS0_C0RE objects that have spectra from previous SDSS surveys that were never 
visually inspected. To guarantee consistency with prior observations, we will visually inspect each of those. The vast 
number of quasar spectra in eBOSS require robust spectral classification in an automated fashion. However, even with 
the new redshifting algorithm presented in Section 5.2, 7% of eBOSS quasar targets will need visual inspection. Given 
that these will be the most difficult spectra to characterize, the time investment will not be terribly different from 
that in BOSS. Improvements to the spectral templates are therefore required to improve the quasar classification in 
an automated fashion. Improvements to the LRG and ELG templates should also mitigate catastrophic failures and 
provide better discrimination between the true redshift and other local minima in the distribution of fits. In addition 
to redshifts for clustering measurements, TDSS and SPIDERS will produce a more heterogeneous spectroscopic sample 
than is expected for the eBOSS LRG, ELG, and quasar clustering samples. Thus new templates will likely be needed to 
meet the requirements for automated redshift determination for the LRG, ELG, and quasar samples, and to encompass 
the heterogeneity of the TDSS and SPIDERS data. 

The BOSS spectrold pipeline was based upon PGA templates for quasars and galaxies and physically motivated 
templates for stars. The BOSS quasar templates do not cover the full variation of actual quasars for BOSS, much less the 
TDSS and SPIDERS samples. Unphysical combinations of the BOSS PGA templates can also lead to erroneous redshift 
estimates as shown in Section 5.1. We are exploring a purely archetype-based alternative where, by construction, every 
template is physically meaningful, thus eliminating the false redshift estimates from unphysical combinations of PGA- 
based templates. We also expect to introduce more linear combinations of templates to allow classification of such as 
binary stars or AGN with significant galaxy light. 

Any new developments in the spectrold component of the data reduction pipeline must be backward compatible 
with SDSS-I, SDSS-II, and BOSS. This requirement is necessary to ensure that redshifts from all samples including 
objects observed in SDSS and BOSS can be determined in a consistent fashion for clustering measurements. 

The initial updates of spectral templates will be modeled on the SEQUELS data to more accurately represent the 
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Table 8 

Expected volume density of each target class, presented in units 10“"^/r®Mpc“^. Entries highlighted in bold font denote the fraction of 
the sample that satisfies the high-level requirement for the redshift distriloution of the sample. Because Lyo quasars will be used as 
random sightlines to projected neutral hydrogen, we report the surface density of those objects parenthetically. 


Redshift 

CMASS 

eBOSS 

Clustering 

Lya 

Fisher 

Low Density 

High Density 


LRG 

LRG 

Quasars^ 

Quasars'^ 

ELG 

DECam ELG 

DECam ELG 

0.6 < z < 0.7 

1.137 

0.810 

0.119 


1.412 

0.183 

0.205 

0.7 < z < 0.8 

0.170 

0.678 

0.130 


2.165 

1.908 

2.068 

0.8 < z < 0.9 

0.010 

0.350 

0.154 


1.654 

2.673 

3.034 

0.9 < z < 1.0 

0.001 

0.097 

0.171 


0.624 

1.135 

1.605 

1.0 < z < 1.1 



0.163 


0.218 

0.373 

0.568 

1.1 < z < 1.2 



0.170 


0.081 

0.159 

0.241 

1.2 < z < 1.4 



0.175 





1.4 < z < 1.6 



0.166 





1.6 < z < 1.8 



0.151 





1.8 < z < 2.0 



0.137 





2.0 < z < 2.1 



0.122 





2.1 < z < 2.2 



0.093 

0.069 (2.75 deg-2) 




2.2 < z < 2.4 




0.063 (5.02 deg-2) 




2.4 < z < 2.6 




0.042 (3.35 deg-2) 




2.6 < z < 2.8 




0.027 (2.16 deg-2) 




2.8 < z < 3.0 




0.018 G-39 deg-2) 




3.0 < z < 3.5 




0.009 G-69 deg-2) 




Observed Surface Area 

7000 deg^ 

7000 deg^ 

7500 deg^ 

7500 deg^ 

1500 deg^ 

1400 deg^ 

1100 deg^ 

Uniform Surface Area 

- 

6120 deg2 

6750 deg2 

- 

- 

- 

- 


^The clustering quasars are uniformly selected and include BAL and non-BAL quasars previously observed in SDSS-I, -II, and -III. 

*^The Lya forest quasars are selected by various means and only include non-BAL quasars that will be observed in eBOSS, hence the 
reduction in surface density relative to clustering quasars over 2.1 < z < 2.2. 

true target spectral variations. Redshifts for LRGs, ELGs, and quasars will be determined by visual inspection and 
used as input in the generation of templates. Subsequent updates of these templates will be possible as more spectra 
are taken (especially of rarer objects) and any problems in the templates are uncovered. We expect to present new 
templates in a future publication similar in philosophy to Bolton et al. (2012). 

7. COSMOLOGICAL PROJECTIONS 

Following the methodology developed in Font-Ribera et al. (2014b), we forecast cosmological constraints for the 
eBOSS program. We begin by using the expected density, redshift distribution, volume, and bias for each tracer to 
independently predict the constraints on the matter power spectrum at the relevant redshift. From this, we predict the 
sensitivity of eBOSS to measure the cosmic distance scale through BAO and to measure the rate of structure growth 
through RSD. We conclude by presenting the power of the combined sample to constrain the properties of dark energy, 
gravity, neutrino masses, and inflation. 

7.1. Number Density of Spectroscopic Sample 

The density and redshift distributions used in projections are taken from the parent target samples described 
in Table 1 and Table 2. Estimates for the number densities of the final spectroscopic sample are corrected for 
incompleteness introduced in the tiling process as explained in Section 4.3. The assumed volume density for each 
spectroscopic sample for the appropriate redshift ranges and the expected survey area is shown in Table 8. In addition 
to reporting the observed area, we report the total area that meets the requirement for 15% uniformity in target 
density. As discussed in Section 3.2, areas with larger fluctuations could be discarded in the cosmology analysis. 
Because we currently have no way of tracking the regions with excess fluctuations, we do not include the loss of 
area due to calibration uncertainty. Instead, zeropoint errors will need to be modeled in the analysis. Likewise, the 
members of the eBOSS collaboration are developing a new bright star mask using the WISE photometry. While the 
final mask has not been established, it is likely that 3-4% of the sky will be removed due to proximity to bright stars 
and suspect photometry. Throughout, a flat AGDM model is assumed with Ho=70 and Dm = 0.31. These predicted 
number densities allow us to predict the sensitivity to Da{z), H{z), and fas from each class of tracer. 

We assume 100% completeness on the clustering quasars observed prior to eBOSS, 95% completeness on the new 
targeted clustering quasars, and 100% completeness on the reobservations of known Lya quasars. For the PTF-selected 
quasars, we assume 92% completeness in fiber assignments and then reduce the density by an additional 10% to account 
for incomplete PTF coverage of the eBOSS footprint. We simply average the effect of incomplete coverage over the 
full area that will be observed. 

The surface density of the LRG target sample is taken to be the more conservative estimate presented in Table 1. 
Because the effective area of the LRG program is reduced by 7% due to masking by targets tiled in the first round, 
we take the average density of targets to be the total number of galaxies assigned a fiber over the reduced footprint. 
We assume a completeness of 87% of the LRG target sample over a final area of ~ 7000 deg^ for the LRG sample 
instead of 7500 deg^. To account for isolated regions with decollided completeness below 85% (see Section 4.3), we 
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include a 5% loss of area in addition to the 8% area that could be lost to targeting non-uniformities when reporting 
the “Uniform Surface Area” for the LRG sample. Finally, because the high redshift tail overlaps with the new LRG 
sample, we include the n{z) distribution for GMASS in the redshift range of interest. In the cosmological projections 
that follow, we combine the two samples for optimal constraints over the redshift range 0.6 < z < 1.0. 

The ELG program is not yet finalized so we include the volume density for three potential selection schemes. We 
do not report the results for the gri + Uri selection as it produces tracers at an efficiency of only 52.5%. We assume 
a completeness of 95% over 1500 deg^ for the Fisher, over 1400 deg^ for the low density DEGam selections, and over 
1100 deg^ for the high density DEGam selection. As it produces a sample with the highest median redshift, we assume 
the high density DEGam-based ELG selection in the projected BAG and RSD measurements and for the cosmological 
projections presented in Section 7.3. 


7.2. BAO and RSD Forecasts 

In order to delineate clustering measurements between samples and ensure the predictions are independent, we use 
only the 0.9 < 2 < 2.2 component of the quasar sample for direct clustering predictions and only the z > 2.1 component 
of the quasar sample for Lya forest predictions. In reality, the QS0_C0RE sample will extend over both redshift regimes 
and will allow for direct clustering measurements at redshifts beyond those used in these projections. 

The Lya quasar sample in eBOSS will complement the BOSS survey, improving the BAO measurement by providing 
deeper spectra of known z > 2.1 quasars and spectra of new z > 2.1 quasars. As explained in Font-Ribera et al. (2014a), 
the auto-correlation Lya forest clustering measurement is shot-noise limited and therefore nearly uncorrelated with 
the Lya forest - quasar cross-correlation measurement. We therefore use the combined precision from these two 
measurements to forecast the final BAO distance measurement. Rather than applying direct Fisher projections, we 
instead scale the expected DR12 results according to the scheme presented in Section 5.3. We do this for H{z) and 
dA^z) but not R because the weighting of the radial modes and transverse modes for determining the optimal distance 
indicator is so different between the cross-correlation and auto-correlation analysis. The measurement of RSD from 
the Lya sample is considerably more challenging for both theoretical and instrumental considerations, and we do not 
project its precision in this document. 

We present the projections for each sample in two year increments, according to the data release schedule presented 
in Section 4.5. For BAO forecasts we assume measurements using modes with k < 0.5/iMpc“^; for RSD forecasts we 
use modes with k < 0.2/iMpc“^. This sequence of projections roughly corresponds to the expected schedule for public 
data releases and the likely publications of clustering measurements. For the LRG, quasar, and Lya forest samples, 
we assume that eBOSS will complete 3000 deg^ in the first release, an additional 1500 deg^ for the second release, 
and an additional 3000 deg^ in the final release. The usable LRG area is somewhat less than the observed area due 
to loss of targets in regions that are masked by higher priority TDSS, SPIDERS, and quasar targets. We report the 
area that is expected to be observed in each case. For the ELG sample, we only present the results expected after the 
second release, in which 300 plates are expected to be observed. The assumed bias for each sample is the same as that 
presented in the beginning of Section 3. 

The time series of predictions on BAO distance precision and RSD growth precision are shown in Table 9. After 
six years, the LRG sample is expected to produce the highest redshift, sub-percent precision distance measurement 
obtained by any means. The quasar sample is expected to produce a 1.8% BAO distance measurement, the first 
percent-level distance measurement achieved in the interval 1 < z < 2. The ELG sample will produce the highest 
precision BAO measurement to date using this tracer of the matter density field. If one were to assume that the 
area that fails to meet 15% uniformity in target selection cannot be used for clustering measurements, the precision 
of the BAO and RSD measurements will be degraded. For the LRG sample, reducing the area to account for low 
completeness sectors and excessive systematic fluctuations results in a 7% increase in the projected errors. For the 
clustering quasar sample, the potential loss of area would result in a 5.5% increase in the projected errors. 

The six-year predictions for the eBOSS distance measurements are shown in Figure 1. It can be seen in the bottom 
panel that SNe la offer high precision measurements of distance out to redshifts z < 0.6 or so, while the BAO probes 
offer constraints comparable to SNe la around z ~ 0.6 but offer much stronger constraints for redshifts beyond. As 

Table 9 

Basic parameters expected for each eBOSS sample, together with predictions for the effective volumes and fractional constraints on BAO 

distance measurements and growth of structure. 


Sample 

Epoch 

Area (deg 

CFH jH 

OD a! Da 

(trIR 

^ fCTR / /^8 

LRG 

year 2 

2790 

0.032 

0.017 

0.012 

0.040 


year 4 

4185 

0.026 

0.015 

0.010 

0.034 


year 6 

6975 

0.021 

0.012 

0.008 

0.026 

ELG (High Density DECam) 

year 4 

1100 

0.047 

0.031 

0.020 

0.038 

Quasar 

year 2 

3000 

0.066 

0.043 

0.028 

0.050 


year 4 

4500 

0.054 

0.036 

0.023 

0.041 


year 6 

7500 

0.042 

0.028 

0.018 

0.032 

BOSS Lyo: Quasars 


10,400 

0.02 

0.025 

- 

- 

BOSS + eBOSS 

year 2 

3000 

0.017 

0.021 

- 

- 

Lyo; Quasars 

year 4 

4500 

0.016 

0.020 

- 

- 


year 6 

7500 

0.014 

0.017 

- 

- 
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emphasized in Kim et al. (2015), BAO and SNe provide distinct cosmological constraints even at the same redshift 
because of the difference between absolute and Hubble flow calibration, making the combination much more powerful 
than either probe in isolation. The combination of SNe from DES with the final data sets from SNLS and SDSS-II will 
provide a far larger sample for reducing both statistical and systematic errors. However, because the leading supernova 
analyses are currently limited by systematic uncertainties (e.g. Sullivan et al. 2011; Betoule et al. 2014), it is difficult 
to forecast the performance of future surveys. 

7.3. Cosmological Constraints 

In what follows, we report the projections for cosmological constraints from the eBOSS BAO distance, RSD, and 
power spectrum measurements from the galaxies and quasars using modes with k < O.HiMpc“^. We report predicted 
cosmological constraints for a six year program using the combined information of all four tracers. We assume that 
clustering measurements can be recovered over the areas which display deviations from the uniformity requirement. 
For these projections, we take the values of the parameters of the fiducial cosmology to be the flat ACDM model with 
parameter values as shown in Table 1 of Font-Ribera et al. (2014b). We allow free parameters to describe the growth 
rate of structure (Ay and a normalization constant Gg), the summed neutrino mass (^TOj/), non-Gaussianity in the 
perturbations of the initial density field (/nl), and a time-evolving equation of state for dark energy (wq, Wa)- We 
assume a baseline of Planck CMB measurements, 5% Hq constraint, and the BAO measurements from the complete 
BOSS galaxy sample. The projected constraints are found in Table 10. A brief interpretation of those predicted 
constraints in terms of dark energy, modified gravity, neutrinos and inflation is presented below. 


Table 10 

Predicted precision from the combination of CMB and large-scale structure measurements. All values correspond to the estimated 1—tr 

uncertainties. 


Parameter 

Constraint from 
CMB 

Constraint from 
BOSS and CMB 

Constraint from 

BOSS, eBOSS, and CMB 

VlAlh? 

0.008 

0.0028 

0.0017 

WQ 

0.52 

0.17 

0.15 

Wa 

1.4 

0.67 

0.48 

7 

30. 

0.13 

0.10 


0.81 eV 

0.29 eV 

0.16 eV 

Us 

0.0045 

0.0026 

0.0022 


7.3.1. Dark Energy 

To demonstrate the power of BAO as an independent probe of the cosmic distance ladder, we report the constraints 
on the equation of state for dark energy using BOSS and eBOSS BAO measurements as the sole source of information 
from cosmic times between the CMB and today. We fix neutrino mass and inflation parameters to highlight the 


Planck+H||+BOSS (w/o Lya), FoM=l 1.3 

Planck+H/BOSS (w/Lya), FoM=24.3 

+eBOSS (w/o Lya), FoM=24.9 

+eBOSS (w/Lya), FoM=73.7 

X 

X 

h 



'i.O -1.5 -1.0 -0.5 -2.0 -1.5 -LO -0.5 0.0 


Figure 12. Current constraints on the DETF model for time-varying dark energy compared to projected constraints from 
eBOSS. We report constraints from the BAO probes, Planck, and Hq from HST observations of SNe la and Cepheid variables 
(Riess et al. 2011). For all measurements, the filled ellipse represents the 68% confidence interval and the open ellipse represents 
the 95% confidence interval. 
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relative constraints on dark energy. We predict constraints on a standard parameterization of dark energy in which 
the equation of state evolves with time as w(a) = wq + Wa{l — a). The prediction for the combined constraints on wq 
and Wa is shown in Figure 12. 


7.3.2. Modified Gravity 

At a given redshift, RSD measurements constrain the product but the extended redshift range of the combined 
eBOSS and BOSS measurements will signihcantly reduce the degeneracy between / and a^. In predicting the power 
of eBOSS RSD measurements to constrain models of modified gravity, we assume a parameterization f(z) = 
where 7 is the growth index (Linder 2005). The value of 7 is approximately 0.55 when predicting the growth rate 
from General Relativity. Measurements with RSD that differ from this value would indicate a model for gravity that 
deviates from GR at cosmological scales. The amplitude of f{z) is normalized according to a free parameter referenced 
as Gg in the assumed model. 

An example of the power of RSD is shown in Figure 2, where the growth rates are compared for a set of cosmological 
models that predict the same background expansion (i.e. the same distance-redshift relation and comoving BAG 
position). In the ACDM model, which assumes that GR correctly describes gravity on all scales, the evolution of fas 
is determined entirely by the expansion history H(z). On the other hand, if GR requires modification on cosmological 
scales, as demonstrated in the other two models plotted, then measurements of /erg over a large redshift interval will 
reveal that inconsistency, thus indicating that cosmic acceleration arises from modihed gravity instead of dark energy. 
The eBOSS measurements will significantly improve the constraints on fas between 0.6 < 2 < 2.0, which will help 
constrain the amplitude of the growth rate at fairly early times where various models nearly converge. More generally, 
the combination of BAO and RSD measurements will enable eBOSS to simultaneously constrain models of dark energy 
and modified gravity (e.g., Song & Percival 2009). A review of the power of eBOSS and other spectroscopic surveys 
to measure the growth of structure in tandem with imaging surveys is presented in Huterer et al. (2015). 

7.3.3. Neutrinos 

A global fit to solar and atmospheric neutrino flavor oscillations implies a difference in the squares of masses 
m 2 — m\ = 7 . 54 /q ’22 X 10“® eV^ and |(m| — m\)l2 + (rug — m 2 )/ 2 | = 2.42/g°g x 10“^ eV^ (Fogli et al. 2012). 
Assuming a normal hierarchy in which the lightest neutrino is massless, one can infer the minimum sum of the masses 
to be greater than 0.05 eV. In the case of a so-called inverted hierarchy, where m 3 = 0, the minimum sum of the 
masses must be greater than 0.1 eV. Both of these are well below the reach of the present terrestrial experiments. 
For example, beta-decay experiments currently constrain the effective electron-neutrino mass to be ;< 2eV (Otten 
& Weinheimer 2008) and even the next generation direct tritium decay experiments such as the Karlsruhe Tritium 
Neutrino experiment (KATRIN; KATRIN collaboration 2001; Wolf & KATRIN Gollaboration 2010) will only have a 
sensitivity to constrain mivfi} < 0.2 eV (90% upper limit). When considering the small mass differences derived from 
the oscillation experiments, the projected results from KATRIN can at best constrain each individual mass eigenstate 
to m < 0.2 eV, or < 0.6 eV (90% confidence). 

The large effective volume covered by eBOSS makes it sensitive to neutrino masses. The signature of neutrinos 
appear as modulations of clustering in the same large-scale structure maps used for BAO and RSD constraints. An 
example of the predicted suppression of power by massive neutrinos is shown in Figure 13. eBOSS will therefore place 
tight new constraints on this fundamental quantity without any changes to the BAO survey design. The relative 
impact of eBOSS, CMB and other cosmology programs to improve neutrino constraints is reviewed in Abazajian et al. 
(2015b). 

The best current cosmological constraints on ^ from large-scale clustering result from the combination of Planck 
CMB, CMB polarization, lensing and BAO. As reported in Planck Collaboration et al. (2015a), these combined probes 
produce an upper limit ^ mi, < 0.23 eV (95% confidence) when assuming zero curvature. With eBOSS, we can provide 
comparable constraints if one combines information derived from the shape of the broad-band power spectrum (for 
neutrino masses) with the distance measurements derived from BAO. Estimates from the clustering of BOSS galaxies 
have already shown the potential in spectroscopic surveys to constrain the neutrino mass (e.g. Zhao et al. 2013; Beutler 
et al. 2014a). 

The projections in Table 10 include an estimated precision on the neutrino mass that is highly degenerate with the 
modihed gravity parameters. If one were to assume a hat, ACDM cosmology described by GR on cosmological scales, 
we predict a 95% upper limit '^mi, < 0.104 eV when combining the results from the BOSS, eBOSS, and CMB. The 
reported constraint assumes that that measurements of large-scale modes can be made with the eBOSS galaxies and 
quasars at all scales larger than kmax = 0.1/iMpc“^. If one were to assume measurements of large-scale modes with 
wavelengths up to kmax = 0.2/iMpc“^, the projections for the 95% upper limit improve to '^m^ < 0.072 eV. The 
statistical power of eBOSS can potentially provide an upper limit on neutrino masses comparable to the minimum 
allowed mass in an inverted hierarchy. 

An alternative way to constrain neutrino masses with optical spectroscopy is to use the information from the ID 
hux power spectrum of Lya forest. Measurements of the ID hux power spectrum using BOSS quasars (Palanque- 
Delabrouille et al. 2013b) yield the tightest constraints to date using any cosmological tracer. When combining the 
BOSS measurement with CMB and BAO measurements, the sum of the neutrino masses is constrained with a 95% 
upper limit ^ m^ < 0.15 eV (Palanque-Delabrouille et al. 2015b), improved to ^ < 0.12 eV (Palanque-Delabrouille 

et al. 2015a). We do not include projections for new constraints from ID power on the neutrino mass from eBOSS. 
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Figure 13. The suppression of power at small angular scales due to free-streaming of massive neutrinos. A normal hierarchy 
for neutrino masses is assumed. 

The deeper spectra obtained by reobserving 1.2 deg“^ known quasars and the discovery of 1.2 deg“^ new quasars 
at z > 3 should allow tighter constraints on the redshift evolution of the ID flux power spectrum, and thus tighter 
constraints on the neutrino mass. 


7.3.4. Inflation 

Inflation is the best candidate that we have for a theory that simultaneously established the initial conditions for 
structure formation while producing a homogeneous, nearly flat Universe. Inflation can explain small super-horizon 
fluctuations which are nearly Gaussian and scale-independent. A review of inflation and the potential of upcoming 
cosmology programs to improve inflationary models is found in Abazajian et al. (2015a). 

The eBOSS survey volume provides sensitivity to the standard inflationary parameters such as the spectral index 
of primordial matter fluctuations (ris), dehned according to P{k) = fc"”. Most inflationary models predict a value of 
ris slightly less than one. One possible departure from standard inflationary scenarios is expected to appear as small 
deviations from Gaussian fluctuations in the initial density field. This signature of non-Gaussianity in the primordial 
fluctuations could also be visible in the large-scale structure maps from eBOSS. 

A form of potential non-Gaussianity that arises naturally in many inflationary models (Wands 2010) is the “local” 
form, with a gravitational potential described by 

(Z* = ((-S +/nL (<('g - , (7) 

where (jig is an auxiliary Gaussian held. In the limit of /nl = 0, one recovers the Gaussian case. It has been shown 
(Dalai et al. 2008) that such non-Gaussianity will generate a deviation from the standard prediction in which galaxies 
are a linear tracer of the underlying dark matter held on large scales. This correction scales inversely as the square of 
the wavenumber in the matter power spectrum, and it thus becomes important on the larges scales - precisely in the 
limit where non-linearities in the dark-matter huctuations are negligible. 

We hnd that primordial non-Gaussianities of the local form can be constrained to a precision = 12 (68% 
conhdence) by the combination of eBOSS LRG, ELG, and quasars. This measurement will be independent of the 
current Planck bispectrum limits of /nl = 0.8±5. (68% conhdence) (Planck 2015 results. XXVII; Planck Gollaboration 
et al. 2015b). While we do not have forecasts for constraints arising from measurements of the galaxy bispectrum, we 
expect these to signihcantly tighten the constraints from eBOSS, provided systematic effects can be brought under 
control (e.g. Ross et al. 2013a). In fact, the hrst analysis of the BOSS galaxy bispectrum (Gil-Marfn et al. 2015a) and 
cosmological interpretation (Gil-Marfn et al. 2015b) were completed in 2014. 
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8. CONCLUSION 

As described throughout this paper, the eBOSS survey design has been thoroughly evaluated and should meet a 
high-level goal of precise clustering measurements using four distinct tracers over four different redshift intervals. A 
sample of 52 deg“^ targets will be observed and will produce more than 35 deg“^ LRGs with confident spectroscopic 
redshifts between 0.6 < z < 1.0. When combined with the z > 0.6 tail of the BOSS CMASS sample, this LRG 
sample has the statistical power to constrain the matter power spectrum at BAO scales for better than a 1% distance 
measurement. Similarly, a uniform sample of quasar candidates selected at a density of 115 deg“^ can be efficiently 
identified. Roughly 30 targets deg“^ have confident spectroscopic classification from SDSS-1, -11, or -111, revealing 
quasars at a density 13 deg“^ between 0.9 < z < 2.2. New observations will produce an average of 56 deg“^ new 
quasars at 0.9 < z < 2.2. The final sample should be sufficient for a 1.8% BAO distance measurement. The range 
of redshifts covered by the LRG and quasar program is only sparsely sampled by supernova surveys, but it covers 
the crucial epoch of transition from decelerating cosmic expansion to accelerating expansion. As with BOSS, the area 
covered by the LRG and quasar clustering samples approaches the limits accessible by ground-based telescopes: we 
will observe one-third of the extra-Galactic sky. 

eBOSS will use a deeper sample of Lya forest observations to improve the BOSS BAO distance measurements by 
a factor of 1.44 and will introduce a new sample of highly-biased ELG targets. Several potential ELG programs have 
been identified that can provide BAO measurements at 2% precision. The final ELG program will be decided when the 
imaging data sets are better understood and when uniformity of the target selection is fully assessed. Beyond BAO, 
the clustering in LRG, ELG, quasars, and the Lya forest provide broad-band power spectra to further explore the 
cosmological model, particularly in the realm of modified gravity and neutrino mass estimates. The primary analysis 
effort within the eBOSS collaboration will be to mitigate systematic errors due to the non-uniformities in the target 
selection and improve modeling of the observed power spectrum at small scales to capitalize on this statistical power. 

The cosmological projections presented in Section 7.3 assume a standard cosmological model. Because of the overlap 
in redshift range between the ELG sample and the LRG and quasar samples, the eBOSS data also enable techniques 
such as those introduced in McDonald & Seljak (2009) and Seljak (2009) to combine multiple tracers and reduce the 
effects of sample variance. Projections for /nl and RSD from eBOSS following the multi-tracer technique are found in 
the work by Zhao et al. (2015). Zhao et al. (2015) also present Eisher forecasts for the dark energy equation-of-state 
and modified gravity based on a principal component analysis parameterization. 

eBOSS will provide unique spectroscopic information that will complement imaging cosmological surveys such as 
DES. Most imaging probes of dark energy rely on photometric redshift estimates which must be calibrated to extremely 
high accuracy to avoid degradation in dark energy inference. This calibration is typically done using spectroscopically 
calibrated redshifts of objects spanning the photometric properties of the sample of interest. However, existing spectro¬ 
scopic surveys are inherently limited by the difficulty of securing redshifts for faint, high redshift galaxies. Spectroscopy 
from eBOSS offers an alternative approach to photometric redshift calibration through cross-correlation techniques 
(Newman 2008; Menard et al. 2013). The observed degree of correlation between a well-measured spectroscopic sample 
and a photometrically-selected sample provides information on the fraction of the imaging sample at that redshift. 
The possibility of cross-correlation as a tool to calibrate photometric redshifts for DES was the primary reason that 
the 500 deg^ eboss2 region was tiled in the first year (see Figure 3). 

Improvements in the automated data reduction pipeline are underway to ensure that we meet the redshift efficiencies 
required to meet the assumed cosmological precision. The introduction of new spectral templates that cover a more 
limited parameter space should significantly improve the ability of the pipeline to differentiate the best fit redshift 
from contaminating interlopers. When complete, the classifications will produce the largest sample of z > 0.6 LRGs to 
date, a high-redshift ELG sample exceeding that of the WiggleZ Dark Energy Survey (Drinkwater et al. 2010; Blake 
et al. 2011a), a sample of quasars exceeding that of all previous observations combined, and an enhanced sample of 
Lya forest spectra. These new samples will allow studies of galaxy evolution and quasar astrophysics well beyond the 
cosmological studies that drove the survey design. 

Several studies of galaxy evolution and quasar astrophysics have already begun with the early release of the SEQUELS 
sample in DR12 and the first eBOSS data taken in Fall 2014. At the redshifts of the LRG and quasar samples, many 
absorption lines are accessible in the eBOSS spectra. Measurements of quasar absorption features in the vicinity 
of BOSS spectroscopic galaxies revealed the Mg II distribution surrounding LRGs at redshifts as low as z = 0.5 
(Zhu et al. 2014). Composite spectra constructed from Lya forest absorbers shed light on circumgalactic regions 
and the intergalactic medium (Fieri et al. 2014). These studies will be extended to higher redshift with the eBOSS 
data. Studies of the correlations in the Lya forest with other systems are underway following early studies in BOSS 
(e.g. Font-Ribera et al. 2012), while others have the shown the potential of eBOSS to offer new probes of large-scale 
structure and BAO (e.g. Fieri 2014). The composite spectrum of ELGs from limited pilot observations already shows 
rich spectral features such as resonant absorption and non-resonant emission. The first science result from the eBOSS 
program resulted from studies of this composite ELG spectrum (Zhu et al. 2015). These features are rare in other 
spectroscopic samples; eBOSS will allow systematic investigation of these features to explore the effect of gas processes 
on galaxy evolution. The overlap in spectroscopic area covered by eBOSS and wide-field far infrared imaging from 
the Herschel Space Telescope (Filbratt et al. 2010) will enable an exploration of the far-IR luminosity function of the 
quasar host galaxies. Finally, the quasar sample will greatly exceed all prior work both in total numbers and in the 
coverage of luminosity-redshift phase space. The sample can be used to enhance BOSS constraints on the luminosity 
function (e.g. Ross et al. 2013b; McGreer et al. 2013), provide halo occupation statistics that can be used to constrain 
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the duty cycle of quasars as a function of halo mass and quasar luminosity, and to explore the redshift and luminosity 
evolution of quasars. When complete, the eBOSS spectroscopic sample will results in a diverse range of findings both 
for cosmology and for galaxy and quasar science. 
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